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Abstract 

In this paper we consider the averaging of local field-theoretic Poisson brackets in the multi- 
dimensional case. As a result, we construct a local Poisson bracket for the regular Whitham 
system in the multidimensional situation. The procedure is based on the procedure of averaging 
of local conservation laws and follows the Dubrovin - Novikov scheme of the bracket averag- 
ing suggested in one-dimensional case. However, the features of the phase space of modulated 
parameters in higher dimensions lead to a different natural class of the averaged brackets in 
comparison with the one-dimensional situation. Here we suggest a direct procedure of construc- 
tion of the bracket for the Whitham system for d > 1 and discuss the conditions of applicability 
of the corresponding scheme. At the end, we discuss canonical forms of the averaged Poisson 
bracket in the multidimensional case. 

1 Introduction. 

We consider here the Whitham averaging method ([35j [36j [37]) associated with slow modulations of 
periodic or quasiperiodic m-phase solutions of nonlinear systems 

^■(^Vt, ¥>*i,..., ¥>*<»,•••) =° > i = l,. <p = ((p 1 ,...,? 71 ) (1.1) 
which are usually represented in the form 



y? l (x, t) = $ l (ki(U) x 1 + ... + k d (U) x d + w(U) t + O , U) (1.2) 

We are going to consider here systems with d spatial variables (x 1 , . . . , x d ) and one time variable 
t. In these notations the functions k 3 (U) and w(U) play the role of the "wave numbers" and 
"frequencies" of m-phase solutions, while the parameters 6$ represent the "initial phase shifts". The 
parameters U = (U 1 , . . . , U ) can be chosen in an arbitrary way, we just assume that they do not 
change under shifts of the initial phases of solutions Q . 

The functions $*(#) satisfy the system 

^($,0;^,^^,...,^^,...) = , z = l,...,n (1.3) 
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and we have to choose for each value of U some function <fr(0,U) as having "zero initial phase 
shift". The corresponding set of m-phase solutions of fll.l[) can be then represented in the form 
(11.21) . For m-phase solutions of f ll.ip we have in this case k p (U) = (/Cp(U), . . . , fc™(U)), u>(U) = 
(a; 1 (U), . . . , u m (XJ)), 6 Q = (9$, . . . , 9™), where U = (U 1 , . . . , U N ) are the parameters of a solution. 
We will require also that all the functions $*(0,U) are 27r-periodic with respect to each 9 a , a = 
l,...,m. 

Consider a set A of functions <f>(0 + O ,U), depending smoothly on the parameters U and 
satisfying system ( 11. 3p for all U. 

In the Whitham approach the parameters U and 0o become slowly varying functions of x and t: 
U = U(X,T), O = 0o(X,T), where x = (x\...,x d ), X = (X\...,X d ), X« = ex\ T = et (e -»■ 0). 

For the construction of the corresponding asymptotic solution the functions U(X, T) must satisfy 
some system of differential equations (the Whitham system). In the simplest case (see |20j), we try 
to find asymptotic solutions 

<A0,X,T,e) = *(*) + 0, X, t\ e k (1.4) 

fc>0 ^ e ' 

with 27r-periodic in functions VtaA satisfying system (11.11) . i.e. 

F'(cp,e(p T ,eip x i,...,ecp X d,...) = , z = l,...,n 
The function S(X, T) = (5' 1 (X, T), . . . , S' m (X, T)) is called the "modulated phase" of solution 

Assume now that the function ^( O )(0, X, T) belongs to the family A of m-phase solutions of (II. ip 
for all X and T. We have then 

* (o) (0,X,T) = $(0 + o (X,r),U(X,T)) (1.5) 

and 

S$(X,T) = co a (V(X,T)) , S XP (X,T) = k«{U(X,T)) 

as follows after the substitution of (jl.4j) into system fll.ip . 

In the simplest case the functions *&(k){9, X, T) are determined from the linear systems 

L} [U)0o] (X,r)^ fc) (0,X,T) = /^(0,X,T) (1.6) 

where L\jj g ](X, T) is a linear operator defined by the linearization of system (11.31) on the solution 
(11.51) . The solubility conditions of systems fll.6p in the space of periodic functions can be written 
as the conditions of orthogonality of the functions i^(0,X,T) to all the "left eigenvectors" (the 
eigenvectors of the adjoint operator) of the operator If^ g JX,T) corresponding to zero eigenvalue. 

We should say, however, that the solubility conditions of systems (11.61) can actually be quite 
complicated in general multi-phase case, since the eigenspaces of the operators L[u,0 o ] anc ^ ^[v e ] 
on the space of 27r-periodic functions can be rather nontrivial in the multi-phase situation. Thus, 
even the dimension of the kernels of £[u,0 o ] and Ltj g ^ can depend in a highly nontrivial way on the 
values of U. In general, the picture arising in the U-space can be rather complicated. As a result, the 
determination of the next corrections from systems (11.61) is impossible in general multiphase situation 
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and the corrections to the main approximation (II. 5p have more complicated and rather nontrivial 
form (PEE]). 

These difficulties do not arise commonly in the single-phase situation (m = 1) where the behavior 
of eigenvectors of L[\j,e \ and Z| U0o j, as a rule, is quite regular. The solubility conditions of system 
(TUD for k = 1 

L} [Ufflb] (X,T)*f 1) (0,X,r) = / ( \)(^,X,T) (1.7) 

with the relations 

k P T = wip , k pX i = hxp , p,l = l,...,d 

define in this case the Whitham system for the single-phase solutions of (11. ip which plays the central 
role in considering the slow modulations. 

For the multi-phase solutions the Whitham system is usually defined by the orthogonality condi- 
tions of the right-hand part of (11.71) to the maximal set of "regular" left eigenvectors corresponding 
to zero eigenvalues which are defined for all values of U and depend smoothly on U. 

As mentioned above, the construction of asymptotic series (11.4D in the multi-phase case is impos- 
sible in general situation (see [6j [TJ [8]). Nevertheless, the Whitham system HI . 13[) and the leading 
term of the expansion (11.41) play the major role in consideration of modulated solutions also in this 
case, representing the main approximation for the corresponding modulated solutions. The correc- 
tions to the main term have in general more nontrivial form than (jl.4p . but they also tend to zero 
in the limit e^O(00|). 

Let us give here just some incomplete list of the classical papers devoted to the foundations of 
the Whitham method [H SI El El d EJ El HU [121 LITl HH UHl [20l [291 ETil EH [351 ESI [37|. We will be 
interested here only in Hamiltonian aspects of the multi-dimensional Whitham method. In general, 
the article represents the technique and methods developed in the paper [25], applied to the multi- 
dimensional case. However, as we will see, the features of the phase space of modulated parameters 
in higher dimensions lead to a different natural class of the averaged brackets in comparison with 
the one-dimensional situation. 

In the remaining part of the Introduction we will give the definition of "regular" Whitham system 
for complete regular family of m-phase solutions which will be used everywhere below. 

Let us use for simplicity the notation A both for the family of the functions <&(0 + 0q, U) and the 
corresponding family of m-phase solutions of system (II. ip . such that we will denote by A both the 
parameter-dependent families of the 27r-periodic in all 9 a functions 3>(0 + Oq,TJ) and <^u,0 o ]( x ) — 
<fr(kg(U) x q + 6q, U). We will assume everywhere below that the family A represents a smooth family 
of m-phase solutions of system (11. ip in the sense discussed above. 

It is generally assumed that the parameters , u a are independent on the family A, such that 
the full family of the m-phase solutions of (II. ip depends on N = (d + l)m + s, (s > 0) parameters 
U v and m initials phase shifts 8q . In this case it is convenient to represent the parameters U in 
the form U = (k 1; . . . , k^, uj, n), where k p represents the wave numbers, uj - the frequencies of the 
m-phase solutions and n = (n 1 , . . . ,n s ) - some additional parameters (if any). 

It is easy to see that the functions & 9a (6 + 6 , ki, . . . , k^, uj, n), a = l,...,m, 
& n i{0 + 9o, ki, . . . , k d , uj, n), I — 1, . . . , s, belong to the kernel of the operator L\ ki kd ^ n 9 i. In the 
regular case it is natural to assume that the set of the functions (<Jv, ^V) represents the maximal 
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linearly independent set of the kernel vectors of the operator L regularly depending on the parame- 
ters (ki, . . . , k^, n )- For the construction of the "regular" Whitham system we have to require the 
following property of regularity and completeness of the family of m-phase solutions of system (II. 1ft : 

Definition 1.1. 

We call a family A a complete regular family of m-phase solutions of system U.l\) if: 

1) The values k p = (k^, . . . , k™) , uf = (uj 1 , . . . , u m ) represent independent parameters on the 
family A, such that the total set of parameters of the m-phase solutions can be represented in the 
form (U, O ) = (ki, . . . , k d , w, n, 6 Q ); 

2) The functions <&0<*(0 + 6 , ki, . . . , k d , u?, n), & n i(0 + O , k l5 . . . , k d , u>, n) are linearly inde- 
pendent and give the maximal linearly independent set among the kernel vectors of the operator 
U- r ki kd w n i ; smoothly depending on the parameters (ki, . . . , k^, u>, n) on the whole set of param- 
eters; 

3) The operator L 1 -^ k d wne ] ^ as exac ^V m + s linearly independent left eigenvectors with zero 
eigenvalue 

ic[g,(0 + o ) = 4!,..,k^,n](^ + ^o) , q = l,...,m + s 

among the vectors smoothly depending on the parameters (ki, . . . , k^, uj, n) on the whole set of pa- 
rameters. 

By definition, we will call the regular Whitham system for a complete regular family of m-phase 
solutions of (II. ip the conditions of orthogonality of the discrepancy fm(0, X, T) to the functions 



p2tt pAtt d m 9 

/•••/ ^ M]i (H0o(X,T))/| 1) (P,T)-- = O , g = l,..., m + s (1.8) 
with the compatibility conditions 

Kt = u xp , k px> = k ixv , a = 1, . . . , m , p, I = 1, . . . , d (1.9) 

System (II ,8p - (II .9p gives md{d + 1)/2 + (m + s) conditions at every X and T for the parameters 
of the zero approximation \&(o)(0, X, T). 

It is well known that the Whitham system does not include the parameters 9q (X, T) and provides 
restrictions only to the parameters U U (K,T) of the zero approximation. Let us prove here a simple 
lemma which confirms this property under the conditions formulated above. 

Lemma 1.1. 

Under the regularity conditions formulated above the orthogonality conditions U.8\) do not contain 
the functions 8q (X, T) and give constraints only to the functions U V (K,T), having the form 

C^{V)U V T - D^{V)U V XV = , q = l,...,m + s (1.10) 
(with some functions Ci q \\J), # P (U) ; v = 1, . . . , JV, p = 1, . . . , d). 
Proof. 
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Let us write down the part of the function f^), which contains the derivatives 0q T (X, T) and 



6g XP (X,T). We have 



dF i 



d(p J t 



(* ( o), .-.) 2a; Q (X,T)^ 



1 \Tf-? flP _ 

I 1 //Virtfl 1/flT 1 



3)8 

' (o)e/ 3 °ot 



dF 



dp 



, (*(o)> •••) *( 0) ^C- 



^ j \ - " ■ - - / - ~ - (O)0<*0 : v 0T 



ft? 

Ur\ r T 



dF 1 
dpi 



(* (0 ), -..) 2A£(X,r)¥ 



0fL - 



tt w r xPx l 

Let us choose the parameters U in the form 



(O)0 a 0P v ox 



U = (fcJ,...,A£V 



n 



We can write then 
/(i)(0,X,T) 

+ 



i} r + ^u),...) + i;A + °o, u) 



aw/ 3 



7 0T 



+ 



(j F l ($(0 + 6> , U), . . . ) + L) A $ J (0 + O , U) 



00 
u oxp 



The total derivatives dF % jduj^ and dF % jdk^ are identically equal to zero on A according to (II. 3p . 
We have then 



2vr /.27T 



(9) 

U(X,T)] i 



i (0 + o (x,r))/g ) (0,x,T) 



(27r)"> 



since k [u(xt)](^ + ^o(X, T)) are left eigenvectors of L with the zero eigenvalue. 

It is not difficult to see also that all the (X.,T) in the arguments of <fr and disappear after 
the integration with respect to 6, so we get the statement of the Lemma. 

Lemma 1.1 is proved. 



Conditions ( II. 8p together with the compatibility conditions 



k 



pT 



IjJ 



XP 



a = 1 m 



p = l,...,d 



give (m + s) + md = m(d + 1) + s restrictions on the functions U(X, T) which is exactly equal to 
the number of the parameters U u . 
Let us call the system 



CP{U)UZ = D\?> P (U)U XP , q 
a — 1 m , 



1, . . . , m + s 



;i.ii) 



K pT 



•'JfP ' Ul — -L, ...,!/(. , P— 1, ...,d 

the evolutionary part of a regular Whitham system for a complete regular family A, while the 
restrictions 



v P x l 



Hxp 



a = 1, . . . , m 



p,l = l,...,d 



;i.i2) 



will be considered as additional constraints for the evolutionary system ( 11. lip . It is easy to see 
that the constraints ( I1.12p are compatible with the evolutionary system ( II. lip in the sense that the 
restrictions (jl,12p are conserved by system (II. lip being imposed at the initial time. 
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In generic case the derivatives U!£ can be expressed in terms of U^ p from system (11. lip and the 
evolutionary part of a regular Whitham system can be written in the form 

U V T = V?(U)U£, (1.13) 

Thus, system f l 1 . 1 3 [) represents a homogeneous quasi-linear system of Hydrodynamic Type con- 
necting the derivatives of the slow modulated parameters. 



2 Lagrangian and Hamiltonian formulations of the Whitham 
method. 

Together with the formulation of Whitham's method the Lagrangian structure of the equations of 
slow modulations was proposed ([251 ESI EZ])- The method of averaging of Lagrangian function 
introduced by Whitham can be formulated in the following way. We assume that the original system 
(II. ip is lagrangian with the local action of the form 



S = J L {if, (ft, ip x , (f tt , (p xt , y> xx , . . . ) d xdt 

such that the functions F l have the form 

i 5S dL d dL d dL 

Let us assume for simplicity that the parameters (k p , a>) = (k^, . . . , k™, u 1 , . . . , u m ) give the 
complete set of independent parameters on the family of m-phase solutions (excluding the initial 
phase shifts), such that the number of parameters U u is equal to m(d + 1). 

The linearized operator ^-r kl kd w e i(X, T) in (II. 6p is given now by the distribution 

5 2 S 

L )[ki,...,k d ,u,e ](OiQ') 



where 



S = J J L[9, u a <S> ea , fcf 1 *^,...,^"*^, 



d m e 

(27r)"-' 



and is a self-adjoint operator. 

Throughout the paper we will always understand the integration with respect to 6 as the averaging 
procedure. For this reason, all the integrals over d m 9 will be defined with the factor l/(27r) m . In 
particular, we will also assume that the variation derivatives of the type SS/S(p l (6) are defined in the 
way that 







on the space of 27r-periodic in functions. 

We also define here the delta function 5(6 — 6') and its higher derivatives (Wi...^ (0 — 8') on the 
space of 27r-periodic functions by the formula 
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/2tt p2tt d m 9' 

The functions 3>6> q , ol = l,...,m represent both the left and the right eigenfunctions of the 
operator £}r kl kd w o i(X, T), corresponding to the zero eigenvalue. 

Under the assumption that the family of the m-phase solutions A is a complete regular family 
of m-phase solutions of (11. ip we assume that the functions + 6 ,ki, . . . ,k d ,uj) are linearly 

independent and give the maximal linearly independent set among the kernel vectors of the operator 
L*[ ki kd u 0q ] smoothly depending on the parameters (ki, . . . , k^, uj). The regular Whitham system 
is given then by the conditions k^ T = ou XP , k® xl = kf XP , and m conditions of orthogonality of the 
function f(i)(0, X, T) to the functions &q<x(0 + O , k 1; . . . , k d , uj). 

According to the Whitham procedure the Whitham system on the parameters (k 1; . . . , k^, uj) is 
obtained from the condition of extremality of the action 



p p Z7r pZ-k (l rn f) 

E<°>[S] = J J J L(#, S?*,., ^#^,...,5^*^,...) j-^d A XdT (2.1) 

under the conditions = Sj^p, w a = SjL 

The conditions fcp T = a/£ P , = kf XP , and 5£/5S a (X, T) = give a system of mc?(d + l)/2 + m 
equations on the parameters (kp, uj). 

It is not difficult to see that the system given by the variation of the "averaged" action coin- 
cides with the conditions of orthogonality of the function f(i)(0, X, T) to the functions <Ev(0 + 
#0) ki, . . . , ka, uj). Indeed, let us consider the action 

E[s ,,, £l = J L ^(^ + e^T)^ v ^ + e^T)^)^x d r = 

= [S, <p] + e £« [S, <p) + e 2 [S, <p] + . . . 

defined on the functions <p{Q, X, T), 27r-periodic in each 8 a . Taking into account the relation 

^ -r>r.T*(w">™ (2.2) 



6S«(X,T) J J ™ v ' ' ' <^(0,X,T) (2tt) 

and the invariance of the action with respect to the shifts 

S(X,T) S(X,T) + AS 

it is easy to see that 

^«(f,X,JJ - = U (2.d) 



o ./o 



5^(0, X,T) (2?r) 



^^^X,T) - (2.4) 



<5^(X,T) 7 ---y ™ V ' ' 7 M»,X,T) (2tt) 

etc. 
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Substituting the functions <p(0, X, T) in the form ip(0, X, T) = <fr(6 + , S x i, • • • , S X d, &t) m the 
relations above, we can see that we have to include now the additional dependence of the functions 
(p(0, X, T) on Sx and St in relation (|2,4p . However, due to the relation 



<V<(9,X,T) 

on the family A, relation (12. 4p will not change in this situation. Taking also into account the equality 

r (1) (o,x,T) = 5{pl{e ^ T) 

we get the required statement. 

Under the assumption of the completeness and regularity of the family A we can see then that 
the averaged action (12.11) defines a lagrangian structure of the regular Whitham system in general 
multiphase case. We should say also that the cases with additional parameters n, as a rule, can be 
also included into the scheme described above with the aid of the Whitham "pseudo-phases" ([37]). 

In this paper we are going to consider the Hamiltonian formulation of the Whitham equations. 
Let us consider then another approach to the construction of the regular Whitham system which is 
connected with the method of averaging of conservation laws. According to further consideration of 
the Hamiltonian structure of the Whitham equations we will assume now that system (II. ip is written 
in an evolutionary form 

<f] = F i (ip, cp x , y> xx , . . . ) (2.5) 
The families of the m-phase solutions of (12. 5p are defined then by solutions of the system 

co a ^ a = F l (ip, <p eh <p 90d , . . . ) (2.6) 

on the space of 27r-periodic in each 8 a functions . 

We will assume that the conservation laws of system (I2.5P have the form 

P? {<p, <p x , <p„, . . .) = Q u x \(<p, v?x, ¥>xx,-..) + ••• + 0${<P, <£x, (2.7) 
such that the values 

F = J P "(^ ¥>*, d d x 

represent translationally invariant conservative quantities for system (12. 5p in the case of the rapidly 
decreasing at infinity functions v ? ( x )- We can also define the conservation laws for system (12. 5p in 
the periodic case with fixed periods K\, . . . , Kj. 

l v = TT \ •••/ P" (¥>, ¥>x, ¥>xx, • • • ) d d x 



K x ...K d j j Q 
or in the quasiperiodic case 



r = & wy £■ ■ ■ C p " (v ' v " v - - ■ ' A 
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It is natural also to define the variation derivatives of the functionals l v with respect to the 
variations of <^(x) having the same periodic or quasiperiodic properties as the original functions. 
Easy to see then that the standard Euler - Lagrange expressions for the variation derivatives can be 
used in this case. 

Let us write the functionals I u in the general form 



V = J P V (<P, ¥>x, ¥>*x,...) d d x 



assuming the appropriate definition in the corresponding situations. 

Let us define a quasiperiodic function <£>(x) with fixed quasiperiods (ki, 
on M. d coming from a smooth periodic function <p(0) on the torus T m : 



(2.i 



kj) as function ip(x) 



^(kix 1 + • ■ ■ + \t d x d + O ) ->• <p(x\ 
Let us define the functionals 

-2tt r 2n 



X 



/Z7T /*Z7T 



d m e 

(27r) Tn 



(2.9) 



on the space of 27r-periodic in functions. 
It's not difficult to see that the functions 



5.P 



8^(0) 



(2.10) 



V5(0) = *(0+0 O ,U) 



represent left eigenvectors of the operator LJ, ue , with zero eigenvalues, regularly depending on 
parameters U on a fixed smooth family A. 

Indeed, the operator Ljp^, is defined in this case by the distribution 



We have 



[6,6') = 8)u a 8 0a (0-0') 



<p(0)=<S>(0+0 ,U) 



2tt r 2n 



8J V 



Pa 



d m e 



o 



,„ J 6^(0) V" rw " ' ^'" x ^'■■■'-«* twv-jj (2?r)m 

for any translationally invariant integral of (12. 5p . Taking the variation derivative of this relation with 
respect to <^(0') on A we get the required statement. 
Thus, we can write 



m+s 



V 



N 



(2.11) 



9=1 



with some smooth functions c v A\J) on a complete regular family A. 



For our consideration of the regular Whitham system we will need a sufficient number of the first 
integrals (12.81) . such that the values of the functionals J u on A represent the full set of parameters 
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U u = J v \\. Thus, we will require here the presence of N = m(d + 1) + s independent integrals I v ', 
v — 1, . . . , m(d + 1) + s. Besides that, we will require that the maximal linearly independent subset 
of the functions (12.101) gives a complete set of linearly independent left eigenvectors of the operator 
L l j\xj O ] with zero eigenvalues among the vectors regularly depending on the parameters U on the 
family A. 

Coming back to the definition of a complete regular family of m-phase solutions of system (12.51) 
we can see that in the case of a complete regular family A the number of linearly independent vectors 
(I2.10p on A is always finite. More precisely, if N = m(d + 1) + s is the number of parameters of 
m-phase solutions of (12. 51) (excluding the initial phase shifts) then for a complete regular family of 
m-phase solutions we require the presence of exactly m + s = N — md left eigenvectors i^^(d + #o) 
with zero eigenvalues, regularly depending on parameters (in accordance with the number of the 
vectors ^V, ^V)- Thus, according to Definition 1.1, we assume here that the number of linearly 
independent vectors defined by formula (12.101) is exactly equal to m + s = N — md for a complete 
regular family A. 

We should note that the conditions on the variation derivatives of J v formulated above do not 
contradict to the condition that the values J u [y = 1, . . . , N) can be chosen as parameters U u on 
the family of m-phase solutions. Indeed, the definition of J u (12.91) explicitly includes the additional 
md functions hp, which provide the necessary functional independence of the values of J u on A. In 
other words, we can use the Euler - Lagrange expressions for the variation derivatives of l v only 
on subspaces with fixed quasiperiods (ki, . . . , k^). The variation of the quasiperiods gives linearly 
growing variations which do not allow to use the Euler - Lagrange expressions. 

Moreover, under the assumptions formulated above, we can show that the condition of the com- 
pleteness of the variation derivatives (I2.10p of the functionals J v in the space of regular left eigenvec- 
tors of the operator Lj ue , with zero eigenvalues follows in fact from the condition that the values 
U u = J u \a can be chosen as the full set of parameters (excluding the initial phase shifts) on the 
family A. 

Let us make the agreement that we will always assume here that the Jacobian of the coordinate 
transformation 

(k ?) u>,n) (U\...,U N ) 

is different from zero on A whenever we say that the values U u (k q , lo, n) represent a complete set of 
parameters on A (excluding the initial phase shifts). 

Under the conditions formulated above let us prove here the following proposition: 

Proposition 2.1. 

Let A be a complete regular family of m-phase solutions of system Ii2.5\) . Let the values (U 1 , . . . , U ) 
of the functionals (J 1 , . . . , J ) Ii2.9\) give a complete set of parameters on A excluding the initial phase 
shifts. Then: 

1) The set of the vectors 

{$^(0 + O , k 9 , oj, n) , $ n i(0 + O , k 9 , w, n) , a = l,...,m, I = 1, . . . , s} 
is linearly independent on A; 

2) The variation derivatives Qui(^ + @o)> given by \2.1(ft) , generate the full space of the regular 
left eigenvectors of the operator L 1 ^ e ^ with zero eigenvalues on the family A. 
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Proof. 

Indeed, we require that the rows given by the derivatives 

dU 1 dU N \ fdU 1 dU N 
a^'"''a^y ' \dn f, "' , ~dn r 

are linearly independent on A. Using the expressions 

dU u r2* f 2n d mQ 



dn l 



r-2-K r2-n rl m R 

J j * A o, u)_ , , = i 



a = 1 m 



on A, we get that the set {<&uja, is linearly independent on A and the number of linearly 

independent variation derivatives (I2.10p is not less than m + s. 

We obtain then that the variation derivatives (12.101) generate in this case a space of regular left 
eigenvectors of the operator LJp 9 , with zero eigenvalues of dimension (m + s). 

Proposition 2.1 is proved. 

Let us prove here the following lemma, which we will need in further considerations. 

Lemma 2.1. 

Let the values U u of the junctionals J u on a complete regular family of m-phase solutions A be 
functionally independent and give a complete set of parameters (excluding initial phase shifts) on A, 
such that we have = k^U 1 , . . . , U N ). Then the functionals /c° (J 1 , . . . , J N ) have zero variation 
derivatives on A. 

Proof. 

As we have seen, the conditions of the Lemma imply the existence of md independent relations 



6r 



i/=i 

on A. 



V(0) 



= , r = l,...,md (2.12) 

ip(9) = *(0+0o,U) 



For the corresponding coordinates U u on A this implies the relations 

N dm 

X)A;(u)dtr = EE^S)( U X( U ) 

v=l p=l 0=1 

for some matrix /i^s(U). 



ucu uo uunoiLici unc niaiiiA A*(^p)(U) as a mdxmd matrix giving a linear transformation p, : M. md — > 
M md between the spaces with bases parametrized by the pairs (/3p) and the index r respectively. Since 
U v provide coordinates on A the matrix /i^(U) has the full rank and, therefore, invertible. We can 
then write 

md N 



T=l U=l 
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The assertion of the Lemma follows then from (j2.12p . 

Lemma 2.1 is proved. 

Let us consider now the system 

(P")t = (Q ul )xi + ■■■ + (Q ud )x* , u = l,...,N = m(d+l) + s (2.13) 

on the space of the functions U(X), where (...) denotes the averaging operation on A defined by 
the formula 



i>2tt p2tt , % d m Q 



,K ' (2tt) 



Let us prove here the following lemma about the connection between systems (12.13)) and (11. lip . 
Lemma 2.2. 

Let the values U v of the junctionals J v on a complete regular family of m-phase solutions A be 
functionally independent and give a complete set of parameters on A excluding the initial phase shifts. 
Then on the space of functions U(X) satisfying the system of constraints 

^r(U(X)) = 3fcr(U(X)) (2U] 
dX l dXP 1 ' 1 

system \2.1§) is equivalent to the evolutionary part U.ll\) of the regular Whitham system. 
Proof. 

Let us introduce the functions 

nr ( ^W,...) = *"- ) , lu-..J d >0 (2.15) 

Using the expression for the evolution of the densities P u (<p, ec£> x > • • • ) we can write the following 
identities 



P?(<p, e^ x , . . . ) = ^'" +h ^ ih - ld \v, ^x, • • • ) {F*{<p, e^ x , • • • )) hxK .. ldXd = 

h,—,ld 

= e^^e^x,-) + ••• + eQ^e^x,...) (2.16) 
To calculate the values e(Q ul ) x i + . . . + t(Q ud )x d let us put now 

(p\e,X) = $*^?21 + 0,U(X)J (2.17) 

where S XP = A£(U(X)). 

The operators ed/dX p acting on the functions (12.171) can be naturally represented as a sum of 
kpd/d9 a and the terms proportional to e. So, any expression f((p,eipx.,---) on the submanifold 
( 12 . 1 Tj) can be naturally represented in the form 

Knap*,...) = £y/R,[*,u] 

l>0 
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where /m[$,U] are smooth functions of ($,4^,$^,...) and (U, Ux, Uxx, • • • ), polynomial in 
the derivatives (Ux, Uxx, • • • ), and having degree I in terms of the total number of derivations of 
U w.r.t. X. Note also that the functions $ appear in fm with the phase shift S(X)/e according 
to ( 12.171) . The common phase shift is not important for the integration with respect to 0, so let us 
assume below that the phase shift S(X)/e is omitted after taking all the differentiations with respect 
to X. 

According to (12. 16)) and ( 12. 6 j) we can write 

2w r2w < ^ d m e 



t(Q ul )xi + ... + e(Q vd ) xd = e J J (Q& m + ... + Qxa 

_ , / / ( Yiv{h...l d ) rpi T[ v(l l ...l d ) w i 

~ J '"J Q ^ \ *[°] r hXK..l d X*[l] "T" U i[l] r hXK..l d X«[U\J / 27r V 
h,...,l d 

= « A f D fas- 4 " - tf-tf ,. -„. i + 

a (7 



(2vr)" 



/■■■/ e (*?- *?* • • • k ~* ■ ■ ■ k > ,i ^ F k + 

t/ 7 I \ 

'1, •■■!»<£ 



1 d ~<1 



- Wyi ii iIn i ti sf/j/j fi,_i)yi .; J yd rnl w y<i ii i[ni l <i ^0/3 z, ..rz^-iiyd rni " 



UiXi...^-!)^ 

d m e 

The last two terms in the above expression represent the integral of the value 



+ U ii i[0] ^ePhX^...l d X d [l] + U u j[l] ^9f>hXK..l d X d \ 



ill — >'d 

and are equal to zero. 

It is not difficult to see also that for arbitrary dependence of parameters U of T, the derivative 
of the average (P v ) w.r.t. T can be written as: 



= jfT E ^ 

Now, we can write the relations (P u )t — (Q^x 1 — ■ ■ ■ — {Q ud )x d = as 
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f- - /^(c^x)]^) [^(e,u(x)) - ift(*,x)] + 

v */ 



a U 



+ l fc iT u x i) 2^ n no] ... k d ...k d ^^l^l^ ^f 

l\,...,l d 

4- (ifi ytHh-ld) , ill t?h 1,-yf t?t-^ &i 

^.y'' 2- 1A i[0] «- fc l ^ W...^...^...^-! / (27T) 



h,...,ld 

where the values Cjru(x)](^) are §i ven (12.101) . 

Consider the convolution (in v) of the above expression with the values dk^/dU v . The expressions 

§§(U(X)) cg ( x„(e) 

are identically equal to zero according to Lemma 2.1. 
From the other hand we have 



■ ■ ■ I yj h kl 1 . . . k 1 1 

Jo Jo , , 



*;! . . . 



X 



X 



n 



v(h...l d ) 



d m e 

(2vr) m 



5 



0A£ 



i- 



°/3 °p 



(2.18) 

<p(0)=*(0,U) 

since the variations of the functions $ are insignificant for the values of according to Lemma 2.1. 

We get then that conditions (12.131) imply the relations k^ T = which are the second part of 
system ( 11. lip . 

Now the conditions 



277 



2tt 



Cg (X )]( ) [^(0,U(X)) 



(2^)""' 



express the conditions of orthogonality of the vectors (I2.10p to the function — $ T 



which 



coincides exactly with the right-hand part of equation (11. 7\i in our case. Since the linear span of the 
vectors (I2.10p coincides with the linear span of the complete set of the regular left eigenvectors of 
the operator L l ^ v 0o i(X, T) with zero eigenvalues, we get that system (12.13P is equivalent to system 
(jl. lip under the constraints (I2.14p . 

Lemma 2.2 is proved. 

Let us note here that it follows from Lemma 2.2 that systems ( 12.131) . obtained from different sets 
of conservation laws are equivalent to each other on the submanifold given by constraints ( 12. 14ft . In 



14 



other words, if system (I2.5P has additional conservation laws of the form (12. 8p then their averaging 
gives relations following from system (12 . 13[) under the additional constraints ( I2.14p . 

The Hamiltonian properties of systems ( 12 .131) and more general systems ( 11. 13ft are very well de- 
veloped in the one-dimensional situation. The general theory of the one- dimensional systems ( 11.131) . 
which are Hamiltonian with respect to local Poisson brackets of Hydrodynamic type (Dubrovin - 
Novikov brackets) was constructed by B.A. Dubrovin and S.P. Novikov. Let us give here a brief de- 
scription of the Dubrovin - Novikov Hamiltonian structures and of the properties of the corresponding 
systems C I 1 . 1 3 f) . 

The Dubrovin - Novikov bracket on the space of fields (U 1 (X), . . . , U (X)) has the form 



{U V {X),U»{Y)} = g^(V)5'(X -Y) + b^(V)Ul5(X - Y) , u,n = l,...,N (2.19) 
The Hamiltonian operator corresponding to (12. 19[) can be written in the form 

> = <T(U)^- + b^(V)Ul (2.20) 

As was shown by B.A. Dubrovin and S.P. Novikov ([91 [TTJ [12]), the expression ( 12 . 1 9[) with non- 
degenerate tensor g u ^{\J) defines a Poisson bracket on the space of fields U(X) if and only if: 

1) Tensor g u ^{\J) gives a symmetric flat pseudo-Riemannian metric with upper indexes on the 
space of parameters (U 1 , . . . , U N ); 

2) The values 



r ^ 7 - V 

where g uX (XJ) gx^iXJ) = 5" represent the Christoffel symbols for the corresponding metric g v/1 (XJ). 



As follows from the statements above, every Dubrovin - Novikov bracket with non-degenerate 
tensor g u ^{\J) can be written in the canonical form ([9| [TTl IT2] ) : 

{n v (X),n^(Y)} = e u 5^5'{X-Y) , e v = ±1 

after the transition to the flat coordinates n u = n"(U) for the metric g Vfl (XJ). 
The functionals 

/+oo 
n v {X)dX 
-oo 

represent the annihilators of the Dubrovin - Novikov bracket while the functional 



/+oo 1 N 
-^(nn\X)dX 
-°° / u=l 



represents the momentum functional for the bracket (I2.19p . 

The Hamiltonian functions in the theory of brackets (12.191) are represented by the functionals of 
Hydrodynamic Type, i.e. 
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/+oo 
h{U)dX (2.21) 
-oo 

The bracket (12. 19ft has also two other important forms on the space of U(X). One of them is the 
"Liouville" form ([9j [TTJ [12]) having the form 



{U»(X),U"(Y)} = ( 7 ^(U) + y-'(U)) 5'(X - Y) + -7T7-T U x 5(X-Y) 



for some functions 7^(U). 

The "Liouville" form of the Dubrovin - Novikov bracket is called also the physical form and 
corresponds to the case when the integrals of coordinates U u 



/+oo 
U"(X)dX 
-oo 



commute with each other. 

Another important form of the Dubrovin - Novikov bracket is the diagonal form. It corresponds 
to the case when the coordinates U v represent the diagonal coordinates for the metric g UII (XJ) and 
the tensor g Ufl (XJ) in (12.19ft has a diagonal form. This form of the Dubrovin - Novikov bracket is 
closely connected with the integration theory of systems of Hydrodynamic Type 

W T = V;{\5)U X (2.22) 
which can be written in the diagonal form 

U% = V V {\J)U V X (2.23) 

(no summation) and are Hamiltonian with respect to the bracket (12.191) . 

It was conjectured by S.P. Novikov that all the systems of Hydrodynamic Type having the form 
(I2.23P and Hamiltonian with respect to any bracket (12.191) are integrable. This conjecture was proved 
by S.P. Tsarev in [31] where the method (the "generalized hodograph method") of integration of these 
systems was suggested. However, the method of Tsarev proved to be applicable to a wider class of 
diagonalizable systems of hydrodynamic type which was called by Tsarev "semi-Hamiltonian" . As 
it turned out later in the class of "semi-Hamiltonian systems" fall also the systems Hamiltonian 
with respect to generalizations of the Dubrovin - Novikov bracket - the weakly nonlocal Mokhov 
- Ferapontov bracket ([2"6]) and the Ferapontov bracket ([131 EH])- Various aspects of the weakly 
nonlocal brackets of Hydrodynamic Type are discussed in [261 EEil (El HSl [32J [23] . 

Let us describe also the procedure for constructing the Dubrovin - Novikov bracket in the one- 
dimensional case when the original system ( 12. 5ft is Hamiltonian with respect to a local field-theoretic 
bracket 

Wixwm = Y^ B %)(^^---) 5(k) ( x -y) ( 2 - 24 ) 

fc>0 

with the local Hamiltonian of the form 

H = I P H {cp,cp x ,...)dx (2.25) 
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which was suggested by B.A. Dubrovin and S.P. Novikov in [HI [HJ [T2] . 

Method of B.A. Dubrovin and S.P. Novikov is based on the existence of N (equal to the number 
of parameters U u of the family of m-phase solutions of (I2.5P ) local integrals 

I v = J P u (c Pl cp x ,...)dx (2.26) 
which commute with Hamiltonian (I2.25f) and with each other 

{I v ,H} = , = (2.27) 

and can be described as follows: 

We calculate the pairwise Poisson brackets of the densities P v having the form 

{P»(x),P»{y)} = Y, A 7{^^...)^ k \x-y) (2.28) 

fc>0 

where 

A?(<p,<p m ,...) = d x Q^(cp,cp x ,...) (2.29) 

according to (I2.27p . 

The corresponding Dubrovin-Novikov bracket on the space of functions U(X) has the form: 

{U"(X),U»(Y)} = (A7)(XJ)5'(X-Y) + -^LJ-IP X 6(X-Y) (2.30) 

Let us remind that we assume that the parameters U v coincide with the values of the functionals 
I" defined on the corresponding quasiperiodic solutions of the family A: 

ir = (r(x)) 

The Whitham system 

(P V ) T = (Q v ) x , u = l,...,N 
is Hamiltonian with respect to the Dubrovin - Novikov bracket (I2.30P with the Hamiltonian 

/+oo 
(P H ) (U(X)) dX (2.31) 
-oo 

The proof of the Jacobi identity for the bracket (12.301) was suggested in [22] under certain as- 
sumptions about the family of m-phase solutions of (12. 5p . Besides that, it was shown in [21] that the 
Dubrovin - Novikov procedure is compatible with the procedure of averaging of local Lagrangian func- 
tions when carrying out of both the procedures is possible. Let us note also that the generalization 
of the Dubrovin - Novikov procedure for the weakly nonlocal case was proposed in [23J . 

In paper [33] all the local brackets (I2.19p for the Whitham equations for KdV, NLS, and SG 
equations were found. Besides that, in papers [21 E] the hierarchies of the weakly nonlocal Hamiltonian 
structures for the Whitham systems for KdV were represented. 

The most detailed discussion and justification of the Dubrovin - Novikov procedure separately 
for the single-phase and the multiphase cases can be found in [22]. In particular, it was shown that 
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the justification of the procedure is in fact insensitive to the appearance of "resonances" which can 
arise in the multi-phase case, which is the basis for its widespread use in the multiphase situation. 

In this paper we are going to consider multi-dimensional evolution systems ( 12.5ft which are Hamil- 
tonian with respect to a local field-theoretic Poisson bracket 

Mx) , ^( y )} = J2 b I,.mM **>•••) -y 1 )--- sM ( xd - y d ) ( 2 - 32 ) 

h,—,ld 

(hi ■ ■ ■ j Id > 0); with a local Hamiltonian of the form 

H = J P H {ff, tfh* ¥>xx, • • • ) d d x (2.33) 

As we said already, we are going to consider complete regular families A of m-phase solutions 
of system (I2.5P satisfying system (12. 6p with some values of (fcp(U), u a {U)). We assume here that 
the solutions of the family A are parametrized by m(d + 1) + s independent parameters (h®, u a , n l ), 
a = 1, . . . , m, p = 1, . . . , d, I = 1, . . . , s excluding initial phase shifts Q%, and we have a set of 
m(d+ 1) + s independent first integrals I v of the form ( 12. 8 p such that their values on A can be chosen 
as the coordinate system {U u } on A (excluding initial phase shifts). We will assume also, that the 
integrals I u commute with each other and with the Hamiltonian H 

{I u , = , {I v , H} = (2.34) 

according to bracket ( I2.32p . 

Let us note that according to ( I2.34p the flows 

= SFfaifr,...) = Mx),f} (2.35) 

generated by the functionals I v according to bracket ( I2.32p commute with the initial flow (12. 5ft . The 
flows ( 12.351) leave invariant the full families of m-phase solutions of ( 12 .51) as well as the values U v = I v 
of the functionals l v on them. For a complete regular family of m-phase solutions with independent 
parameters (kp, . . . , k™) it's not difficult to show that the flows ( I2.35P generate linear (in time) shifts 
of the phases 8q with some constant frequencies u au (XJ), such that 

S^(<&, fcf 1 *^,...^"*^,...) = u a »(XJ)% a (0,V) (2.36) 

According to Lemma 2.1 we have also that the functionals /c~(I) should generate the zero flows 
on the corresponding family of m-phase solutions of f l2.5f) . For U u coinciding with the values of I v 
on A we get then the relations 

dk a (XJ) 

-^J-u Pv (V) = , a,P = l,...,m, p = l,...,d (2.37) 

By analogy with the one- dimensional case we can try to construct an analogue of the Dubrovin - 
Novikov procedure for d > 1. Let us represent the pairwise Poisson brackets of the densities P"(x), 
P M (y) in the form 
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{p-(x) , ^(y)} = £ A Z.iM **.•••) - • • • - 2/) 
(h,...,i d > 0). 

According to relations (I2.34p we can also write here the relations 

4£o(v>, ¥>*,...) = dxiQ^faip*,...) + ... + ^rt,^,...) 

for some functions (Q^ 1 , ■ ■ ■ , Q Ufld ). 

In the full analogy with the Dubrovin - Novikov procedure we can define the expressions 

{C/"(X) , [/"(Y)} = (^io...o)( X ) " yl ) <K x2 - V 2 ) • • • <5(X d - F d ) + ... 

+ (4£ 01 )(X) - F 1 ) <5(X 2 - F 2 ) . . . - Y d ) + (2.38) 

+ (((T^x 1 + • • • + (Q^ d )x*) SiX 1 -Y 1 ) ... 5(X d - Y d ) 

which gives a skew-symmetric (contravariant) form on the space of functions U(X). 

The theory of the Poisson brackets having the form ( I2.38P in d > 1 dimensions was considered in 
[TOt [27J [28] . As was shown in (TOj [27J [28] , the restrictions on the form of f l2.38|) in d > 1 dimensions are 
much stronger than in the one-dimensional case and in general should be considered by the methods 
of the theory of integrable systems. 

However, the procedure described above does not give in general a Poisson bracket on the space 
of fields U(X) in the multi-dimensional situation, since the expression ( I2.38P does not satisfy the 
Jacobi identity for d > 1 in general case. Fortunately, we don't need a Poisson bracket on the full 
space of fields U(X) for the Hamiltonian formulation of the Whitham approach since the regular 
Whitham system is defined on the submanifold in the space of {U(X)} given by constraints (I2.14p . 
As we will show in the next chapters, expressions ( I2.38P define indeed a Poisson structure after the 
restriction on the submanifold mentioned above which gives a Hamiltonian structure for the regular 
Whitham system for d > 1. 

For the description of the Hamiltonian structure of the regular Whitham system for d > 1 it 
is convenient to introduce a coordinate system on the submanifold corresponding to the Whitham 
solutions. It is most natural to consider then the functions S a (X.), a = 1, . . . , m such that S XP (X.) = 
fep(X) as a part of a coordinate system on the submanifold given by constraints f)2.14p . However, 
the X-derivatives of the functions S a (X.) give just md independent parameters in the space U(X) 
at every given X. To get the full system of coordinates on the submanifold defined by constraints 
(I2.14p we need additional m + s parameters at every point X. It will be convenient for us here to 
take arbitrary m + s parameters ?7 7 (X), 7 = 1, . . . , m+s from the set {U(X)} which are functionally 
independent with the set {k® (U)}. 

Let us note that we put here 7 = 1, . . . ,m + s without any loss of generality. As a necessary 
condition of the functional independence of the set {k® (U), U 1 , . . . , JJ m+s } we have, in particular, 
that the variation derivatives (I2.10p of the corresponding functionals (J 1 , . . . , J m+S ) give a linearly 
independent set on the family A. As a consequence, we can claim that the regular left eigen- vectors 
k\™(0 + 0q) of the operator LJ^^ can be also expressed as linear combinations of the vectors 

Qui (0 + Oq), 7 = 1, . . . , m + s, such that we have 
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m+s 

k% (0 + ) = £ dtp (U) Cg] (0 + Oo) (2.39) 

7=1 

with some functions d^\\J) on A. 

Representing the frequencies u a (XJ) as functions of the new parameters 

u a = co a (S x , U\...,U m+s ) 

we can write the second part of the evolution system ( 11. lip together with the constraints ( j 1 . X 2 1) in 
the form 

5° = u/*(S x , U l ,...,U m+s ) (2.40) 
The remaining part of the regular Whitham system can be written as the relations 

Ul = (Q 7 V + ••• + (Q ld )x« , 7 = 1,...,™+* (2-41) 

where (Q 7P ) = (Q 7P )(S X , U 1 , . . . ,U m+s ). 

As follows from our considerations above the choice of the functionally independent parameters 
£7 7 , 7 = 1, . . . ,m + s (and the corresponding choice of the functionals I 1 ) is unessential for the 
construction. 

We will show below, that the regular Whitham system given by equations (I2.40p - ( 12.41 j) is 
Hamiltonian with respect to the Poisson bracket given by the relations 

{S a (X}, SP(Y)} = , {S Q (X), U\Y)} = u a ~< (S x , U l {X),...,U m+s {X)) 5(X - Y) 

a, (3 = 1, . . . , m, 7 = 1, . . . , m + s, and relations (I2.38P for the skew-symmetric form {?7 7 (X) , U P (Y)} 
restricted to the subset 7,p = l,...,m + s, where 

(AZ.1...0) = (AZ.1...0) (Sx,U\...,U m+s ) , (Q lpp ) = (Q^ p ) (S*,U\...,U m+s ) 
The Hamiltonian function for the regular Whitham system is also local in this case and is equal 

to 

H = J (P H ) (S x , ^(X), . . . , U m+s (X.)) d d X 

Let us note here that although system ( I2.40p - (I2.4ip as well as the Hamiltonian structure de- 
scribed above can be formally defined with the aid of the m + s additional integrals J 7 , the presence 
of the full set {I u , v = 1, . . . , m(d + 1) + s} giving the full set of parameters U on the family A is 
important for the justification of the construction of the Hamiltonian structure in our scheme. The 
rest of equations (I2.13P gives in this situation additional conservation laws for the regular Whitham 
system. Let us note also, that the choice of the functionals {J 7 , 7 = 1, . . . , m + s} is also not impor- 
tant for the construction of the Hamiltonian structure of the regular Whitham system as well as for 
the Whitham system itself^ 

1 We mean here that the Hamiltonian structures obtained with the aid of different subsets {/ 7 , 7 — 1, . . . , rn + s} 
transform into each other after the corresponding change of coordinates. 
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In the next chapter we start the consideration of the averaging procedure for the bracket (12.321) 
giving the Hamiltonian structure for the regular Whitham system (I2.40p - (12 .41 p . 



3 The Dirac restriction and the averaging of the Poisson 
brackets. 

Let us start with the definition of a regular Hamiltonian family of m-phase solutions of system (12. 5p 
and of a complete Hamiltonian set of the functionals (12. 8p . 

Definition 3.1. 

We call family A of m-phase solutions of system Ii2.5\) a regular Hamiltonian family if: 

1) It represents a complete regular family of m-phase solutions of Ii2.5\) in the sense of Definition 



2) The corresponding bracket ^2. 32\) has on A constant number of annihilators N 1 , . . . , N s with 
linearly independent variation derivatives 5N l /5(p l (x.) which coincides with the number of independent 
annihilators in the neighborhood of A. 

Let us say that according to the generalized Darboux theorem we can identify the number of 
the variation derivatives d~N l /S(p l (x) on A with the number of linearly independent quasiperiodic 
solutions vf\x) of the equation 



where v\ (x) have the same quasiperiods as the corresponding functions y?(x) on A. 
Definition 3.2. 

We call a set (J 1 , . . . , I N ) of commuting functionals Ii2.8\) a complete Hamiltonian set on a regular 
Hamiltonian family A of m-phase solutions of system Ii2.5\) if: 

1) The restriction of the functionals (I 1 , . . . ,1) on the quasiperiodic solutions of the family A 
gives a complete set of parameters (U 1 , . . . , U ) on this family; 

2) The Hamiltonian flows generated by (J 1 , . . . , I ) generate on A linear phase shifts of 6 with 
frequencies uj v (XJ), such that 



3) The linear space generated by the variation derivatives 8I V /5(p l (x) on A contains the variation 
derivatives of all the annihilators N q of the bracket Ii2.32\) . such that 



for some smooth functions 7^(U) on the family A. 

Let us note that it follows from Definition 3.2 that if a complete Hamiltonian set of integrals 
(J 1 , . . . , J ) exists for a regular Hamiltonian family A then the number of the additional parameters 
(n 1 , . . . ,n s ) discussed above is equal to the number of annihilators of the bracket (I2.32p . Indeed, 



1.1: 








rk | ^(U)! | 



= m 




N 
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according to Definitions 3.1 and 3.2, the number of the functional l v having linearly independent 
variation derivatives on A exactly equals to m + s, where s is the number of annihilators of the 
bracket (I2.32p . The total number of independent parameters U v on A is then equal to m(d + 1) + s 
due to the wave vectors k", a = 1, . . . , m, p — 1, . . . , d which implies the above assertion. 

As follows from the condition (2) of Definition 3.2 and from the invariance of the functionals 
N l and I u with respect to the flows (I2.35p . the values 7^,(U) can always be chosen independent on 
the initial phase shifts on the family A. The values 81" /5(^*(x)|a are linearly dependent on A, so 
it's natural to choose a complete linearly independent subsystem. Remembering that the variation 
derivatives of J u (I2.10p are linear combinations of the regular left eigenvectors ^[^{0 + do), we can 
write 

= E <( U ) K m WU) a? + Ob) (3.1) 

¥J (x)=*(k p (U)xP+0 o ,U) q=1 



Sip* 



for some smooth functions n q (TJ). The functions n q (XJ) are then uniquely determined on A and we 
have rk \ \n (U)|| = s by Definition 3.1. 

Let us say that the full linearly independent subsystem of the regular left eigen-vectors of the 
operator L 1 ^ e , with zero eigen- values can be given also by the variation derivatives of the functionals 
I 1 , . . . , I m+S on A, which give exactly m + s additional parameters U 1 , . . . , U m+S to the parameters 
kp. Relations ( I2.39P give the connection between the eigenvectors i\fh(0 + 6 ) and Cf[uj(0 + ), 

7 = 1, ... ,771 + S. 

Easy to see then that for a complete Hamiltonian set of integrals J 1 , . . . , I N we have also the 
relations 

rk ||w a7 (U)|| = m, a = l,...,m, 7 = 1, . . . ,m + s (3.2) 
for the frequencies corresponding to the functionals J 1 , . . . , F 



Tm+s 



The construction which we are going to consider is in fact closely related with the construction 
of the Dirac restriction of a Poisson bracket on a submanifold. Let us describe here briefly the Dirac 
procedure. Using the terminology of finite-dimensional spaces we can say that the Dirac restriction of 
a Poisson bracket on a submanifold Af k C A4 n is associated with a special choice of coordinates in the 
neighborhood of the submanifold M k . The coordinates in the neighborhood of the submanifold M k 
are divided into the "coordinates on the submanifold" (U 1 , . . . , U k ) and the constraints (g l , . . . , g n ~ k ) 
which define the submanifold Af k . It is assumed that the submanifold Af k is defined by the conditions 

g*(x) = , i = l,...,n — k 

while the functions ?7 1 (x), . . . , U h (x) 011 JW" play the role of a coordinate system on M k after the 
restriction to this submanifold. 

If the Hamiltonian flows generated by the functions U^{x) leave the submanifold M k invariant, 
i.e. we have 

{^(x),^(x)} = for g(x) = 

then the pairwise Poisson brackets of the functions (x) define a Poisson tensor in the coordinates 
(U l ,...,U k ) after restriction on Af k which is called the Dirac restriction of the Poisson bracket 
{...,...} on the submanifold M k C M n . 
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In general, according to the procedure of Dirac, if we have some constraints g*(x) which define 
the submanifold M k and some functions W (x) which give a coordinate system on M k we must find 
k linear combinations fi\ (U) ^ s (x) at every point of Af k , such that for the functions 

E?(x) = L?(x) + Pt{U)g'{x) , j = l,...,k 

we have the relations 

{^'(x),<f(x)} = 

for g(x) = 0. 

The functions £/ J '(x) take the same values that [P(x) at the points of Af k and we can define the 
Dirac bracket {...,... }d on Af k by the formula 

{U\U>} D = {^(x),^(x)}|^(U) 
The functions Pi (XT) are determined from the linear systems 

{^(x),^(x)}U*^"(U) + {<f(x),[P(x)}i^ = , z = l,...,n-fc 
and we can also write 

{U\W} D = {^(x),tP(x)}k* - /3:(U){/(x),^(x)}U & ^(U) 
for the Dirac bracket on Af k . 

Let us now describe the procedure, which will be considered in our situation. We consider now 
system ( 12. 5ft which is Hamiltonian with respect to some local bracket f 1 2 . 3 2 f) with a local Hamiltonian 
function of the form (I2.33[) . We first introduce the extended space of fields 

<p(x) -> <p{0,x) 

where the functions ip(0, x) are 27r-periodic with respect to each 6 a , and define the extended Poisson 
bracket: 

M0,x),^(0',y)} = J2 <,..^,(v,fffx,..O* Pl V 1 -^---^H^-^)^-»0 (3-3) 
Let us then make the replacement x — y X = ex and introduce the Poisson bracket 



{<p\e, X) , tpf(#, Y)}= " h+ "' +ld B l,-,h)^ e ^x, • • • ) S^XX'-Y 1 ) . . . S^(X d -Y d ) 6(0-6') 

h,...,ld 

(3.4) 

on the space of fields (p(&, X). 

Let us describe now the submanifold K, in the space of the functions <p(0, X) which we are going 
to consider. 
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First, we will assume that the functions cp(0,X) G K. represent functions from the family A of 
the m-phase solutions of (12.51) with some parameters U(X) at every X. 
Second, we impose the relations 

» = «) , a = 1 ,..., m , p , i = 1 ,..., rf ,3,) 

for the functions U(X) parametrizing the solutions <p(0, X) G A from the submanifold /C. 

More precisely, let us choose for simplicity the boundary conditions in the form XJ(X , 0, . . . , 0) — > 
Uo, X 1 — > — oo, such that kf (Uo) = 0u and define the functions S a (X) on K by the formula: 

r xl 

S a (X) = / k?(X 1 ',0,...,0)dX v + (3-6) 

J — oo 

rX 2 rX d 

+ / k%(X\X 2 ',0,...,0)dX 2/ + ... + / k«(X\...,X d -\X d ')dX d ' 
Jo Jo 

We define then the functions <p(0, X) G /C by the formula 

ip\e, x) = & {^^- + e, s x , ^(x), . . . , u m+s (x)^j (3.7) 

where the functions fc^(X) = Sj£ p (X), and the additional parameters (U 1 , . . . , t/ m+s ), introduced in 
the previous chapter, play now the role of parameters on the family A. 

The functions {f2 Q (X), f/ 1 (X), . . . , U m+s {X.)} play the role of a coordinate system on the sub- 
manifold /C. 

We have now to introduce the analogous coordinates in the vicinity of the submanifold /C. 



Let us introduce the functionals J U (X) on the functions ip(0,X) by the formula 

d m 

(27T 



r(X) = /.../ P^e^eVxx,-..) 77T^7 . i/ = l,-..,JV (3.8) 



and consider their values on the functions of the family /C. 
We can write on K: 

J"(X) = V V {X) +^e z J ( ^(X) , z/ = l,...,iV (3.9) 

where J^(X) - are polynomials in the derivatives Ux, Uxx, • • • with coefficients depending on U 
and have grading degree I in terms of total number of derivations with respect to X. 

The higher terms in (I3.9p are not uniquely defined on /C due to relations (I3.5p . Let us assume 
here that the terms J^(X) are chosen in some definite way in every order I > 1. The corresponding 
choice will affect the definition of the functionals U(X) in the vicinity of /C in the higher orders in e 
(I > 1). As we will see, this choice will not be important in our further considerations. 

The transformation (13.91) can then be inverted as a formal series in e, such that we can write 



2 The case when the values kf = are absent on the space of parameters U requires just a simple modification of 
the definition of S(X) which we for simplicity do not consider here. 
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i>i 



(3.10) 



on the functions of the submanifold K. In formula f l3.10|) the functions UYn are functions of J, Jx, 
Jxx, • • • , polynomial in the derivatives Jx, Jxx, • • • , and having degree I in terms of the number of 
derivations w.r.t. X. 

We can define now the functionals U(X) on the whole functional space using the definition of 
the functionals J(X) and relations (I3.10p . 

Let us put now the same boundary conditions for the functionals fc"(U) as before on the whole 
functional space. We can then consider also the functionals S^fU^X), given by (13. 6p . as the function- 
als defined in the vicinity of the submanifold /C using the corresponding definition of the functionals 
U(X). On the submanifold /C we will naturally have the relations 

5« P [U](X) = A£(U(X)) 

However, outside the submanifold /C these relations are in general not true. 
Let us introduce also the constraints g % {6, X) defining the submanifold /C by the conditions 
g l (0, X) = 0, and numbered by the values of 6 and X: 



g\6,X) = ^(0,X) - ( S[U[J J ](X) + 0, S«[XJ[J}}, ^[J](X), f/ m+s [J](X)) (3.11) 

The constraints g l (0, X) are functionals on the whole extended space of fields ip l (0, X) by virtue 
of the corresponding definition of the functionals J"(X). 

The constraints ( 13.111) are not independent since the following relations hold identically for the 
"gradients" 5g i (0,X.)/5(p j (0' ,Y) on the submanifold /C: 



r-2n 



2tt 



5S a (Z) 



6<pi(0,X) 



k W^Y) 



K. 



d m e 

(2tt)" 



d d X = 



a 



, m 



(3.12) 



r-2n 



2tt 



V(0,x) 



K 



d rn e 



d d X = 



7 = 1, . . .,m + s 



(3.13) 



Nevertheless, it will be convenient here not to choose an independent subsystem from (13.111) and 
keep the system of constraints in the form (13. lip keeping in mind the existence of identities (I3.12p - 

USE}. 

Thus, we consider now the values of the functionals [S^X), {7 1 (X), . . . , U m+s {X.), g l {6, X)] as a 
coordinate system in the neighborhood of /C with the relations ( I3.12p - ( 13.131) . The values of the 
functionals [^"(X), f/ 1 (X), . . . , U m+s {X)] will be considered as a coordinate system on K. 

For the procedure of the averaging of bracket ( 12.321) we will need to find the pairwise brackets on 
/C of the functionals, introduced above, according to the bracket ( 13.41) . 



The pairwise Poisson brackets of the functionals J^(X), J M (Y) have the form 
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{J"(X) , J"(Y)} 

"2tt />2tt 



;>27r ;>27r rl m f) 

. , Jo Jo l Z7r J 



where 



AZ (<P,ecp x ,...) = e a 11 Q^ 1 ^x,..) + ...+e5 xi g^,^x,..) (3.14) 
The Poisson brackets of the fields <p l (6, X) with the functionals J M (Y) can be written as 

ill— j'd 

(3-15) 

with some smooth functions CT ; s(<£, e </?x, • • • )• 
We also have in this case 

q o %M0,X),e^ x (0,X),...) = ^(^Xj.e^fl.X),...) (3.16) 

by virtue of (12351) . 

For any function of slow variable g(Y) we can write 

|V(0,X), f q(Y) J^(Y) d d Y\ = £ e l ^-^C^ ld) (cp(e,X),ecp^e,X),...) q hX i. JdX a 
< J J h,...,i d 

(3.17) 

The leading term in expression (I3.17P on /C has the form 



V \0,X), / q(Y)J^Y)d d Y 



K[0] 



^)(*(— + 0.X).... ) ,/<X) 



g(X) 5*" ( $ f^^ + 0,X , . 



\ e 

where S lfJ, (ip, <£> x , . . . ) is the flow (I2.35P generated by the functional J M . According to f !2.36|) we can 
write then: 

X),y g(Y)J^(Y)^Y} = w^(X) (^ + 0,U(X)) g(X) + 0(e) (3.18) 
Similarly, for any smooth function Q(0,X), 27r-periodic in each 8 a , we can write on the basis of 
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JJ..jJ Q (m +e ^ W 0,X),J"(Y)}| K = (3.19) 

p2ir p2ir rl m fj 

= u<*»(X) jf y Q(0,X) $* a (0,U(X)) *(X-Y) + 0(e) 

In view of relations ( 13. 9 p - ( I3.10P for the functionals J M (Y), U^(Y) we can also write 

|^(9, X),y 9 (Y)f/"(Y)d d r| = w°"(X) *"„ ^?<5> +9,U(X)J g (X) + 0( £ ) (3.20) 

y 2 ' . ^ g (520 + e, x) {,-(9, x), u>m}\c = 

/■2-7T p2tt fl rn f) 

= uT»(X) ... Q(6,X) <^(0,U(X)) — — 5(X - Y) + 0(e) 
</o Jo \ Z7T ) 



= O(e) (3.22) 

K. 



In this case, by virtue of (I2.37P we have 

p*(0,X),y g(Y) A; Q (U(Y)) d<V 

y...y Q^-Li + 0,Xj &(d,X) , k"(V(Y))}\ K = 0(e) (3.23) 

for fixed values of coordinates [S(Z), U\Z), . . . , f/ m+s (Z)]. 

Let us prove now some lemmas about structure of the Poisson brackets on the submanifold K, 
which we will need in the further consideration. 

Lemma 3.1. 

Let the values U u of the functionals I u on a complete regular family A of m-phase solutions of 
112.5]) be functionally independent and give a complete set of parameters on A, excluding the initial 
phases. Then for the Poisson brackets of the functionals k^ (U(X)) and J M (Y) on fC we have the 
following relations: 

{^(U(X)),J^(Y)}k = e [c^(U(X)) 5(X-Y)} XP + 0(e 2 ) (3.24) 



Proof. 

The conditions of Lemma 3.1 coincide with the conditions of Lemmas 2.1 and 2.2. Consider the 
Hamiltonian flow generated by the functional J q(Y) J M (Y) d d Y according to bracket (13.41) with a 
compactly supported function q{Y) of the slow variable Y = ey. 

According to (13.171) we can write at the "points" of the submanifold /C: 

<p\ = q(X) uA(U(X)) (^ + 0,U(X)) + e V \ q] ^S. + 0, x) + 0(e 2 ) (3.25) 
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with some (27r-periodic in each 9 a ) functions 77? q A0, X). 

Let us introduce the functional A£(J(X)) = k^J^X), . . . , J N (X)) using the same functions 
fcp(U). According to (13. 9p - (13.101) we can write 

{^(U(X)),J^(Y)}k = {^(J(X)),J^(Y)}k + 0(e 2 ) 

Consider the evolution of the functionals fcp(J(X)) according to the flow generated by 
J q(Y) J M (Y) d d Y on the submanifold /C. Using relations (I3.25P we can write 

dk a 

WW) = ^ (J(x)) jt(x) = 



^ (J(X)) g(X) uA(U(X)) x 

c2-k p2tt 



x /■'■/*£ n ." ( " " (* + ». U(X)) , e A. (?2> + 9 , U(X)) , . . . ) x 



x e 



h+ - +ld -W777T- ■ ■ ■ -W777T- + 0, U(X 



dX ll i "' <9X<^ e V e ' V (2tt) 



+ e ^j(U( X ))E (?(X)^(U(X))) X9 x 

9=1 



/ 2?r . . f* £ l q Uf^ (* (0, U(X)) , ^ ^ (0, U(X)) $ e7d (0, U(X)) , . . . 

Jo Jo i i 



x 



h,...,ld 



+ e^(U(X))x 



/ 2Tr . . f J2 n: (h ~ ld) (# (0, u(x)) , fc? 1 <&<m (0, u(x)) , . . . , kf ^ (o, u(x)) , . . . ) x 



x 



x k?...k< ...k?...k<rii al a} ad af (e,X)^L + (e 2 ) 

h ...e a i...e l d {27r) m 

It's not difficult to see that the first part of the above expression contains the integrals over of 
the expressions P^(0, X) and is equal to zero. It is easy to see also after integration by parts that 
the third part of the above expression represents the value 



e 



2ir p2n f)fca fj m f) 



and is equal to zero according to Lemma 2.1. 
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Thus, we can write in the main order the expressions for the evolution of the functionals (J(X)) 
on the submanifold K, in the form 



^( J ( x )) = e off-, (u(x)) e (g(x) ^(u(x))) xq y . . . y Yl l « u tf ld) 

9=1 ^ ^ h,...,ld 



X 



9Po a i...e 'i ...e a i...e ...e a i...g l d {2ir) m 

Using relations (" 12 . 18[) we get then 

*£(J(X)) = e[g(X)^(U(X))] XP + 0(e 2 ) 

i.e. 

{r(J(X)),^(Y)}|^ = ew^(U(X))^(X-Y) + e^5(X-Y) + 0(e 2 ) 
which implies (I3.24p by virtue of relations 03.9p - (I3.10p . 

Lemma 3.1 is proved. 

Easy to see that according to Lemma 3.1 and relations (13. 9p - (I3.10P we can write also 

{A£(U(X)), IFOOMz = e [w^(U(X)) S(X- Y)] XP + 0(e 2 ) (3.26) 
on the submanifold /C. 

Let us formulate now some corollaries following from Lemma 3.1. 

1) Under the conditions of Lemma 3.1, we have also 

{A£(X) , kf (Y)}\ K = 0(e 2 ) (3.27) 

for the functionals k p (U(X)). 

Indeed, by virtue of (I2.37P we have 

{A£(X) , kf (Y)}\ K = {A£(X),Lr*(Y)}|jc^(U(Y)) = eu^(X)k^(X)6 XP (X-Y) + 
+ eu°*{X) (k? m (X)) 5(X-Y) + e (co^(X)) XP k^(X)S(X-Y) + 0(e 2 ) = 0(e 2 ) 

2) Using definition (13. 6 j) of the functionals 5 a (X) we can write on /C under the conditions of 
Lemma 3.1: 

{S a (X), J»(Y)}\ K = 6 ^(U(X)) 6(X-Y) + 0(e 2 ) (3.28) 
{S Q (X), ^(Y)}|x: = ew a "(U(X)) 5(X-Y) + 0(e 2 ) (3.29) 
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for the functionals S a (X). 



0(e 2 ) 



{S a (X),S^(Y)}\ !C = 0(e 2 ) 



(3.30) 



Let us note now that in the proof of Lemma 3.1, we have not used the fact that the functionals 
J^(Y) belong to our special set of functionals (13.81) and used only the fact that the flow generated by 
the functional 1^ leaves invariant the family A generating linear shifts of 6q with constant frequencies 
u afl (XJ). We can therefore formulate here the following lemma: 

Lemma 3.1'. 

Let the values U u of the functionals l v on a complete regular family A of m-phase solutions of 
system A2.5\) be functionally independent and give a complete set of parameters on A, excluding the 
initial phases. Let the flow generated by the functional 



P(<P, V* 



d d x 



(3.31) 



leave invariant the family A generating linear shifts of 6$ with constant frequencies cD a (U). 
Consider the functionals 



J(X) 



2n r 2* ^ jmg 



o Jo {^) 



(3.32) 



Then for the Poisson brackets of the functionals k a (XJ(X)) and J(Y) on K we have the relation: 
{^(U(X)),J(Y)}k = e[^(V(X))5(X-Y)] XP + 0(e 2 ) 



Proof of Lemma 3.1' completely repeats the proof of Lemma 3.1. 
Lemma 3.2. 

Let the values U u of the functionals l v on a complete regular family A of m-phase solutions of 
system A2.5\) be functionally independent and give a complete set of parameters on A, excluding the 
initial phases. Then for the constraints g i {0, X) imposed by A3.11\) and smooth compactly supported 
function q(X) as well as smooth 2n-periodic in each 6 a function Q(0,X) we have the following 
relations on the submanifold K: 



2tt p2n 



Q 



g l (6,X), I q{Y).P{Y)d d Y 
S(X 



K 



d,X) {g\0,X),J»(Y)} 



0(e) 

oT9 
(2-n) m 



0(e) 



(3.33) 
(3.34) 



K 



Proof. 

Indeed, by A3. 18j) . f)3.19p . and Lemma 3.1 we have 



g\0,X), / q(Y) J"(Y) d d Y 



KM 



SfX) 



+ 0,U(X) u^{X)q{X)- 
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srxi 



+ 0, U(X) <^ S Q (X) , / g(Y) J"(Y) cfV 



/C[Q] 



/ '"q(^ + 0,x') {^(0,x),^(y)} (/ '"" 



Jo 



(27T) 



/C[0] 

p2tt p2tt tj m f) 

- U °*(X) J J Q(0,X) ^(0,U(X)) ^ <5(X — Y) 
Q(0,X) $» a (0,U(X)) {S Q (X), ^(Y)}!^ 



o Jo 



(27T) 



Lemma 3.2 is proved. 



Similarly to the previous case, we can formulate also the following lemma: 
Lemma 3.2'. 

Let the values U u of the junctionals I u on a complete regular family A of m-phase solutions of 
system A2.5\) be functionally independent and give a complete set of parameters on A, excluding the 
initial phases. Let the flow generated by the functional A3. 31\) leave invariant the family A generat- 
ing linear shifts of 9q with constant frequencies u a (XJ). Then for the constraints g l (6,X) and the 
functionals J{X) imposed by A3. 32\) we have the following relations on the submanifold JC: 



2tt r 2n 



Q 



g l (6,X), I q(Y)J(Y)d d Y 

srx 



K 



0,X W(0,X),J(Y) 



0(e) 

d m 6 
(2vrj^ 



0(e) 



(3.35) 
(3.36) 



K 



under the same conditions for the functions q(X) and Q(6, X) 
Using transformations (13.91) - f 1 3 . 1 j) we can also write 



2tt r 2n 



Q 



'0 Jo 
for the functionals f/ M (Y). 



g l (0,X),J q(Y) U"(Y) d d Y 
S(X) 



K 



0,X {g*(0,X),U»(Y)} 



0(e) 

d m e 

(27r) OT 



0(e) 



(3.37) 
(3.38) 



K 



The pairwise Poisson brackets of the constraints g l (6,X), g^(6',Y) on /C can be written in the 
form: 



{g i (0,X),g*(e> t Y)}\ K = 

W(o,x) ,^(e',Y)}\jc - {^(e,x),u\Y)}\ K: % 

'S(X 



(3.39) 



0'.U(Y)] - 



+ 0,u(x) {u v (x),^(e',Y)}\ K + 
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+ %„ (?21 + 0,U(X)) {U V {X),U\Y)}\ K & ux (^p- + 0',V(Y)^ 
- \ M0,X), S\Y)}\ K (^ + 0\U(Y)) - 

1 % a (^- + e,v(x)) {s a (x),^(o',Y)}\ K + 



e \ e 



+ \ ^ (^ + «,U(X)) {S a (X),U\Y)}\ K % x (^p- + 0',V(Y)^j + 

+ -¥ Tr J^ + e,u(^)){u v (X),s^Y)}\ K %J^- + e',iJ(Y)) + 



e 2 



?21 + 0,u(x)) {s a (x),s^Y)}\ K {^p- + e'MY)\ 



Let us note that we assume that the parameters U", U x represent here a full set of parameters U 
(z/, A = 1, . . . , N) on A. The choice of the parameters U u is in fact not important due to the invariance 
of the corresponding expressions with respect to the substitutions U v = U U (XJ) on A. So, for the 
sake of brevity, we put that the set (U 1 , . . . , U N ) represents the set of parameters (k^, U 1 , . . . , JJ m+s ) 
or the set of averaged densities of the functionals I v introduced above. 

For convenience let us introduce the functional 



J W = j J"(Y) q,(Y) d d Y 



for any smooth compactly supported vector-function q(Y) = (gi(Y), . . . , q^f(Y)). According to 
Lemma 3.2 we can write the relations {g l (6,X) , </[qi}|/c = 0(e) on the submanifold /C. More 
precisely, the first non- vanishing term of the Poisson bracket of g l (0, X) and Jm on K. can be written 

as 



{g\e,x),j h] }\ m = Mm.J^km - $J^ + fl,u(x)j{ff(x), MMm- 

- *U + 0,U(X)) {^(X) , J [q ]}k[2] (3.40) 

Let us remind here that the indexes [1], [2] mean as before the terms of the corresponding graded 
expansion on /C having degree 1 and 2 respectively. 

As we mentioned already, the procedure of the averaging of a Poisson bracket is connected to 
some extend with the Dirac procedure of the restriction of a Poisson bracket onto a submanifold. The 
more detailed consideration of the connection of the averaging procedure with the Dirac procedure 
can be found in |25j. For the justification of the averaging procedure considered here we will need to 
investigate the solubility of the system 

B^(X)B M (6,X) + 4 ][q] (0,X) = (3.41) 

where 
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sg(x) = ^WTO^W^^ 11 ^)) TO^lfl.upc)),...) x 

x^(X)...^(X)...^(X)...^(X) ^...JL ..JL..-^ (3,2) 
is the Hamiltonian operator f |2.32j) on the family of m-phase solutions of system (12. 5p and 

4n q ](^ + 0> x ) = Wkm ~ ^(^)+e,u(x)){t/"(x), J [q] }k w 

(3.43) 

is a discrepancy connected with the first non- vanishing term of the Poisson bracket of g*(0,X) and 
J[q] on /C. Let us note also that for the sake of brevity we denote again by U v the full set of parameters 
U on A as in the expression (I3.39[) . 

System (13.411) represents at every X a linear system of partial differential equations in with 
periodic coefficients. For us the solubility of system ( I3.4ip on the space of 27r-periodic in all 6 a 
functions will be important. Let us discuss now the properties of system (I3.4ip . 

It is easy to see that the operator -BE(X) = _B^(U(X)) is a differential operator on the torus. 
It can be seen also that for special values of k p (X) the foliation defined by the set of the vector fields 
(k^X), . . . , k^(X)) can define the tori of the lower dimensions T k C T m and even one-dimensional 
tori embedded in T m . 

The operator £?E(U) has in general finite number of "regular" eigenvectors with zero eigenvalues 
defined for all values of the parameters U and smoothly depending on the parameters. However, for 
special values of U the set of eigenvectors with zero eigenvalues can be infinite and is determined, in 
particular, by the dimension of closures of the foliation leaves in T m , defined by the vectors k p (U)Jf] 

Let us define in the space of the parameters U the set Ai , such that for all U G M. the dimensions 
of the closures of the foliation leaves, defined by the set {k p (U)} in T m is equal to m. From the 
condition 



rk \\dk%/dU u \\ = md 

it follows that the set A4 is everywhere dense in the parameter space U and, moreover, has the full 
measure. We note also that the set U represents here the full set of parameters U = (U 1 , . . . , U N ) 
on A excluding the initial phase shifts 6q. 

3 For some special brackets (|2 .32[) the differential part can be absent in the operator _B^(U). The operator B%A\J) 

reduces then to an ultralocal operator acting independently at every point of T m . As a rule, the matrix BliJTJ) is 
non-degenerate in this case. For example, for the ultralocal Poisson bracket 

{V>(x), $(y)} = iS{x-y) 

for the multi-dimensional NLS equation 

iipt = A ip + v\ijj\ 2 ijj 

we have exactly this situation. Easy to see that system p.41[) represents a simple algebraic system in this case and 
is trivially solvable. The multiphase situation is not different here from the single-phase case. However, for arbitrary 
brackets (|2.32j) the operators BlL (U) have more general form described above. 
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In the study of the solubility of ( I3.4ip we must first require the orthogonality of the functions 
^■[i][q](0> X) to the "regular" eigenvectors of _B^(U(X)) with zero eigenvalues. Let us prove here the 
following lemma. 

Lemma 3.3. 

Let the values U u of the junctionals l v on a complete regular family A of m-phase solutions of 
system k2. 5\) be functionally independent and give a complete set of parameters on A, excluding the 
initial phases. Then we have the relations 

-Jf 4u(x)]W4lH(^ X ) (2^T = ' ? = l,-,m + « (3.44) 
for the functions A!fa<M), X) defined by \3.43\ ). 
Proof. 

Let us first prove the following statement: 

The values Qu(x)](^)> T = 1, • - • , + s, are orthogonal (for any X) to the values A l ^^(0,X.), i.e. 

-jf $(x)](fl)4]w( fl ' X )^ = > l = h...,m + s (3.45) 

Indeed, according to ( 13. 13f) the convolution of the values 5U 1 {Z)/5ip l {0^ X) with the values 
{^(0,X) , J[q]}|jc identically vanish on K. According to O - flBTTUD the values 5W(Z)/5ip i (6,X.) 
coincide in the leading order with the values SJ 1 (Z)/5(f l (0, X) on /C. 

Using the explicit expression for the quantities 5J~ 1 (Z)/d~(p l (6, X) on /C: 



5J"(Z) 



^ n 7Gi-i d ) U H^p- + 0, zY . . . ) e h+ -+ l * 5^{Z l -X 1 ) ... 5^\Z d - X° 



V(0,x) 

we can write the corresponding convolution in the form of action of the operator 



on the distributions {<?*(#, Z) , Jj q ]}|/c, which is given in the main order by the action of the operator 



/>2n /*27r 
'0 JO 



v 1 . . rv d . " J 



(2< 



d d d d 



x k" 1 (Z)...k?"(Z) ... k a /(Z)...k ' d (Z) -^...-^ ... " (3.47) 

The leading order of the brackets {g % (0, Z) , Jj q ]} on /C is given by the value {g l {6, Z) , J[ q ]}|jqi], 
defined by formula f)3.40p . After integration by parts we get that the action of the operator (13.471) is 
given by the convolution w.r.t. 6 of the values {g l {0, Z) , J[q]}|/c[i] with the values Q™L z y,(S(Z)/e+0). 
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We know also that the values Q™L z s,(S(Z)/e + 0) are automatically orthogonal to the functions 
$g a (S(Z)/e + 8, Z), so (after the replacement of Z to X) we get relation f !3.45j) . 
Using now relations (I2.39P we get the statement of the Lemma. 

Lemma 3.3. is proved. 

For a regular Hamiltonian family A and a complete Hamiltonian set of integrals (J 1 , . . . , I N ) we 
can also prove the following lemma: 

Lemma 3.4. 

Let the functions Bj^Ad, X) satisfy conditions flff.^ifj . Then the functions £>,,[ q ](0, X) automati- 
cally satisfy the conditions 

... & 0a (O, S x , U 1 (X),...,U m+s (X))B M (e,X) 7 — - = , a = l,...,m (3.48) 
Jo Jo l Z7r J 

Proof. 

Indeed, the implementation of (13.41 ft implies the conditions 

/2-7T p2lT ct m Q 

■■l <$ W] (*)^(X)B iW (0,X)^ = , 7 = l,...,m + , 

which is equivalent to 

«^(U(X)) / .../ ^ a (0,U(X))S i[q] (0,X) -— = 

view the skew-symmetry of .Bfi(X). 

From the property (13.21) of the subset (I , . . . , J m+S ) of a complete Hamiltonian family of com- 
muting functionals we immediately obtain now relations (13. 48 p . 

Lemma 3.4 is proved. 

In the remaining part of the article the solubility of system (I3.4ip will play the basic role for the 
justification of the main results. 

In what follows we consider separately the single-phase (m = 1) and the multiphase (m > 2) 
cases. 

The following lemma can be formulated for the single-phase case m — 1: 
Lemma 3.5. 

Let A be a regular Hamiltonian family of single-phase solutions of A2.5\) and (J 1 , . . . ,1) be a 
complete Hamiltonian set of the first integrals of the form A2.8\) . Then the functions BjwAQ, X) can 
be found from system 11(3.41 ) and can be written in the form 



B iM (9,X) = /3f p (0,U(X))^(X) + ^(^U(X))^, g(1 (X) + ptf(e,U(X))Vh 

(3.49) 

(summation in n = 1, . . . , N, p,q = l,...,d, 7 = 1, . . . , m + s) with smooth dependence on the 
parameters U(X). 
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Proof. 

System (13.411) in the single-phase case is a system of ordinary differential equations in 9 with a 
skew-symmetric operator J3E(X). It is easy to see also that the right-hand side of system (13.411) has 
the form 

-4 ][q] (0,X) = ^*(0,U(X))^(X)+#"(0,U(X))^ x , ^(X)+£^(0,U(X))t# P g,(X) 

with periodic in 9 functions £^ p (#, U(X)), £^(0, U(X)), £^ p (#, U(X)). 

The orthogonality conditions ( 13.441) imply then the orthogonality of all the sets of functions 
f^(0,U(X)), £™(0,U(X)), £^(0, U(X)) to the functions K ( fu ( x)]W> such that system (I3T4B 
can be split into independent inhomogeneous systems: 

flg(X) /3f p (0,U(X)) = r ,p (0,U(X)) (3.50) 
%X)/?£f (0,U(X)) = e w (^U(X)) (3.51) 

^o]( X )^7^ U ( X )) = C P (^ U ( X )) ( 3 - 52 ) 

defining functions (13.491) . 

All the systems (I3.50p - (13.521) are systems of ordinary linear differential equations with periodic 
coefficients and a skew-symmetric operator B^ (X). The zero modes of the operator B^ (X) are given 
by the variation derivatives of annihilators of the bracket (12.321) and are orthogonal to the right-hand 
parts of (I3.50p - (13.521) according to (13.11) and (I3.44p . Eigenfunctions of -B^(X) form a basis in the 

space of 27r-periodic functions (f{9). Besides that, the nonzero eigenvalues of B^JX.) are separated 
from zero in this case. Thus, the 27r-periodic functions ft ' p (9, U(X)), /3^ pq (9, U(X)), /3f; p (0, U(X)) 
can be found from systems (13.501) - (13.521) up to the variation derivatives of the annihilators of the 
bracket (12.321) . If we impose additional conditions of orthogonality of /3f ,p (0, U(X)), /3^ P9 (0, U(X)), 
Pi^ p (9, U(X)) to the variation derivatives of the annihilators of the bracket fl 2 . 3 2 [) on the manifold 
of single-phase solutions we can suggest a unique procedure of construction of these functions. The 
functions ffl ' p (9, U(X)), f3 p a pq (9, U(X)), 0^ p {6, U(X)) then depend smoothly on the parameters 
U(X), which implies the required properties for the functions Bj^(9, X). 

Lemma 3.5 is proved. 
Let us start now the investigation of system (13.41 j) in the general multi-phase case. 
Lemma 3.6. 

Let A be a regular Hamiltonian family of m-phase solutions of system $2. 51) and (J 1 , . . . , I N ) be a 
complete Hamiltonian set of commuting first integrals of $2.5\) having the form ( fjg.gj) . Let for UgM 

v&(0) = (v?l v] (e),...,vU v] (o)) , 1 = 1,. ..,8 

be the complete set of linearly independent eigenvectors of the operator B^(XJ) on the torus with zero 
eigenvalues, smoothly depending on 0. Then 

1 ) The number of the vectors vSi (#) is equal to the number of annihilators of the bracket $2. 32\) on 
the submanifold of m-phase solutions of Ii2. 5\) ; 
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2) The functions ^^^(0, X) are orthogonal to all the vectors v [u(x)](^)> 



v 



/o Jo 



'»fU(X)]W ^[l][q] 



°) 4ur«,l(*.X) (2^p = ° ' (U(X) EM) 



Proof. 

Consider the values of v|^j(#) on any of the leaves of the foliation, defined by the set {k 9 (U)} 
on the torus T m . According to the definition of regular Hamiltonian family of m-phase solutions of 
(12.51) the corresponding functions ^^(kix 1 + • — h kdX d + 6 ) should be the variation derivatives of 
some linear combination of annihilators of the bracket (I2.32p . We have then for a fixed value of 6 : 



(0 



From relation ( 13.1 j) we have then 



5N P 



P =i 



5ip l (x) 



V5=*(kj xJ+0o,U) 



(0 



;k,(U) + 6 ) = J2J2 ^( U > *o) "?( U ) K ?u] M U ) + o) (3.53) 



p=l q=l 



By definition, for U6jV( the leaves of the foliation, defined by the set {k g (U)}, are everywhere 
dense in T m . Since both the left- and the right-hand parts of ( 13.531) are smooth functions on T m , we 
get then that they coincide on T m . We can put then a 1 (U, O ) = acL(U) and write 



s m+s 

V %]W -EE a UU) <(U) k% (6) (3.54) 

p=l q=l 

It is easy to see also that any linear combination of the form (I3.54p gives a regular eigenvector of 
the operator -BE(U(X)) with zero eigenvalue. 

Statement (2) follows then from relation (13.441) in view of the representation ( 13. 54ft . 

Lemma 3.6 is proved. 

We can see, in particular, that the values ^.mrqi^j X) are orthogonal at any X to all the "regular" 
eigen- vectors fSjfX)] W °^ ^ ne °P era tor B^l (X) corresponding to the zero eigen- value view the regular 
dependence of the values Af 1][q] (0, X) on the parameters U(X) = (S x , ^(X), . . . , U m+s (X.)). 

However, despite the presence of Lemma 3.6, study of system (I3.4ip is much more complicated in 
the multiphase case if compared with the single-phase case. Thus, the presence of "resonances" may 
lead to appearance of small eigenvalues of the operator 5E(U(X)) for some values of the parameters 
U(X). As a result, this circumstance may lead to insolubility of system ( I3.4ip in the space of smooth 
periodic (in all 8 a ) functions for the corresponding values of U(X). The solubility of system (13.41 j) 
is thus determined by the properties of the operator B^(U(X.)) for the given values of U(X). 

The set of the "resonant" values of U, as a rule, has measure zero in the full space of parameters. 
Let us prove here the following Theorem which shows that the procedure of the bracket averaging is 
in fact insensitive to the appearance of the resonant values of U and can be used in most multi-phase 
cases, as well as in the single-phase case. 

Theorem 3.1. 
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Let system \2. 5\) be a local Hamiltonian system generated by the functional \2.33\) in the local field- 
theoretic Hamiltonian structure Ii2.32\) . Let A be a regular Hamiltonian family of m-phase solutions 
of Ii2.5\) and (I 1 , . . . , I N ) be a complete Hamiltonian set of commuting integrals A2.8\) for this family. 



Let the parameter space U of the family A have a dense set S C Ai on which system ( 3. 41 ) is 
solvable in the space of smooth 2tt -periodic in each 8 a functions. Then the bracket 

{S a (X), S P {Y)} = , {S a (X),U^{Y)} = w" 1 (S x , U 1 {X) , . . . , U m+S (X)) 5(X — Y) , 

{IP(X),U"(Y)} = (^... )(Sx, U\X),...,U m+s (X)) 5 x i(X-Y) + ...+ 

+ (AZm) (Sx, ^(X), . . . , U m+s (X)) 5 xd (X - Y) + 

+ [(Q 7PP )(S X , U 1 (X),...,U m+s (X))] XP 5(X-Y) , >y,p=l,...,m + 8 (3.55) 
obtained with the aid of the functionals (J 1 , . . . , I N ), satisfies the Jacobi identity. 
Proof. 

As before, for a smooth compactly supported vector- valued function q(X) = (gi(X), . . . , gjy(X)) 
we define the functional 

J [q] = J J V {X) q v {X) d d X 

(summation in v = 1, . . . , N). 

Then, for arbitrary smooth, compactly supported in X and 27i~-periodic in each 8 a functions 
Q(0, X) = (Qi(0, X), . . . , Q n (0, X)) we define the functionals 



Qi(0,X) = gi(e,X)-^(0,U(X))M^(U(X)) r.. /^(^X)^, T (0',U(X)) 

Jo Jo 

where the matrix ikP 7 (U) is the inverse of the matrix 



d m 6' 
(2tt)" 



2tt /•2tt n rj m f) 

^^(0,U)^ 7 (0,U) ' 



o J o 



which is always defined according to the definition of a complete regular family of m-phase solutions 
of system (12. 5p . 

By definition, the functions Qi(0, X) are local functionals of U(X) 

Qi(0,X) = Q t (0,X,U(X)) 

depending also on the arbitrary fixed functions Q(0,X). Everywhere below we will assume that 
Q(0,X) is a functional of this type defined with some function Q(0,X). 

Easy to see that the values of Qi(0,X) with arbitrary Q(0,X) represent for fixed values of the 
functionals U(Z) all possible smooth, compactly supported in X and 27r-periodic in each 8 a functions 
with the only restriction 
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2-7T r2ir r] m 

Qi(0,X) $» o (0,U(X)) 



'o jo (^) m 
For the functionals Qi(0,X) we define the functionals 



VX 



a — 1, ... ,m 



(3.56) 




<f(0,x) g 4 [-LJ. + o,x 



o jo 



d m e 

(2vr) m 



Now, for fixed functions q(X), p(X), and functional Q(0, X) consider the Jacobi identity of the 
form 



{9[Q] , {J[q\ , J[p]}} + { J [p] > {^[Q] > J [q]}} + { J [q] > { J [p] > fl'IQ] } } = 



(3.57) 



Expanding the values of the brackets {J iy (X) , J^(Y)} in the neighborhood of the sub manifold 
JC, we can write: 

{J-(X) , J»(Y)} = {7"(X) , J»(Y)}I K + 



+ E 



/"27T /*27T 

/... / AZ. ld ^(e',X),e^(6',X),...) 
Jo Jo 



5<p k (0,W) 



d m 9' 
(2vr)" 



x 



x g k (0, W) - — — d d W e ll+ - +ld S^iX 1 - F 1 ) ... 5 {ld) {X d - F d ) + 0(g 2 ) 
(27r) m 



where all the values on the submanifold K, are calculated at the same values of the functionals [U(Z)] 
as for the original function tp(8, X). 
Let us introduce by definition 



6g"(0,W) 



(3.58) 



a: 



d,...,(j 



d m 6' 
(2vr)" 



x 



x g,(X)p MihXl .., dXd (X)^X 

Note that the notations 5/5g k (6,W), generally speaking, are not natural in our situation be- 
cause of the dependence of the chosen system of constraints. Nevertheless, the preservation of these 
notations can better clarify the algebraic structure of the further calculations. 

The values defined by (13 . 58[) can be represented in the form of the graded decompositions at 
e — > on the submanifold /C. By virtue of (I3.14p it is easy to conclude that the expansion in e of the 
quantities (I3.58P begins with the first degree in e 



<K J [q] > j [p]} 



5g k (0,W) 



<K</[q] > J [p]} 



0,W 



+ ^ 



<H J [q] > J [P]} 



K[l] 



0,W 



+ 



(3.59) 



K[2] 
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The leading term of (13.591) can be divided into two parts, corresponding to the sets of the functions 

(eQ^V efe ...),...,eQ>,^x,..)) 

and 

and containing the quantities <2v(W)^(W) and Q , i/(W)p MjH/ i(W) as local factors, respectively. 
The values of {Jj q ] , J[ p ]}, are obviously invariant under transformations of the form 

(p(0,X) ->■ <^(0 + A0,X) (3.60) 

From the invariance of the functionals U(X) in such transformations, we can write for the corre- 
sponding increments of constraints (13. Ill) 

Sg k (6,X) = (p k (0 + 50,X) - (p k (0,X) 
As a consequence, we can write on the submanifold K 

ff 27T t HJi^J\p]} 

J Jo" 'Jo 8g*(e,W) 

(3.61) 

a = 1, . . . , m. 

The above relation is satisfied to all orders in e. By the arbitrariness of the functions ^(X), 
relation (13.611) in the leading order can be strengthened. Namely, according to the remark about the 
form of the leading term of (I3.59P we can write for any W 



/S(W) 



+ 0, S w ,U\W),...,U m+s (W] 



d m 9 

(27T) 7 " 



d d W 



f 2n f 2n &{J[<A>J\P]} 
l Q "'J 5g k (6,W) 



K[l) 



(^p- + 0, S w , U\W), . . . , U m+s (W)^j = (3.62) 



a = 1, . . . ,m. 



At the same time we have the relations 
<K J [q]> J [p]} 



K 



5S a (W) 

on the submanifold /C. 

Quite similarly, we have the relation 



0(e) 



<H J [q]> J [p]} 



0(e) 



K 



d d X + 



{#[Q] > J[q\} = jQi(^p- + 0,x) {^(0,X), J [q] } I 

+ j g \e,x)Q i>e J^ + e,x\ i {s»(x), j [q] } ^rf d x + 



+ Jg\9,X)Q hU „(^ + 0,X^ {lT(X),J M } ^d d X 



40 



Qj?3p. + o,x\ ^ 1+ -- + ' d qT,., d )M^x) )eV x(0,x),... 

\ / I- I , 



ll,—,ld 



d m e 

(271-)""' 



Q f j,,i 1 x 1 ...i d x d 



(X) d d X - 



ft |?ffl +9iX WS(X) 



g, + e. x i * 



in 



+ 0, Sx^^X),...,^^^ {S£,(X), J [q] } 
^ + ! S x ,f/ 1 (X),...,[/ m+s (X)) {[F(X),J M } — 



(27r) m 



d d X - 
d d X - 



0, x) $*„ + e,s x ,u 1 (x),...,u 
J g \e,x)Q^n*fQ + e,x\ 1 {5»(x), j w } (2?r 



i r/ m $ 

— X)] -{s*(x), j w } d"x + 



d m 



Let us note that the summation in 7 is made here for 7 = 1, . . . , m + s, while the summation in 
v is made for the full set of U u , v = 1, . . . , N. 

According to relations (I3.56p . we can actually see here that the fourth term in the above expression 
is identically equal to zero. 

According to the form of the constraints, we have again in the neighborhood of /C: 

{9[Qh J [q]} = id[Q], J[q\}\,c + 



; X^ e h+-+ld 

h,—,ld 



^ Jd) (#,X), e¥Jx (0',X),...) d m 9 > 



Qti,hx 1 ...i d x d (X-) d d X 



x 



K 



d m e 



x 9 k (0,W) ^ d d W - 



^ S(X) +0',X) ®i„(^W + 0', S x ,U\X),...,U m+s (X)\ S^x 



{S%p(X) , J [q ]} 



V=(0,W) 



(27T) 

d m f) 

/ X >< (2^^ 



- J Qi(0',X) & w (0', S x ,U\X),...,U m+s (X)) 



d m 8> 5{C/T(X),J [q] } 



(2vr) m 6<p k (0,W) 



d d X x 



imfl 

x ^>W)^^ + 
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d m 9 
(2vr) m 



d d W + 



+ jQ k ,e*(^p- + 0,W^ l - {S a (W), J [qi }\ K x g\e,W) 

+ J Qi,u*(^p- + 0, W^j {U^(W) , J [n] }\ K x g k (0,W)-^d d W + 0(g 2 ) 

provided that all the values on the submanifold K are calculated at the same values of the functional 
[U(Z)]. 

For the quantities { <7[q] , J\q\ } we can introduce by definition 



5g k (6, W) 



(3.63) 



/S(X) 



0',X x 



*Cg w (^',X),6^x(^X),... 



ilH Hd (*l-..*d) 



V=(0, W) 



Q f i,hx 1 ...i d x d 



X) 



a: 



(2vr)" 



- Q 



/S(X) 



+ 0', xj $1„ f + 0', S x , ^(X), . . . , f/ m +*(X)) 



cT0' {^ P (X),J [q] } 



(27r) m 5(p k (0,W) 



d d X- 



- J Qi(o\x) (e f , Sx^x),...,^^)) 



cf"0' 5{C/T(X),J [q] } 



(27r) m 5ip k (e,W) 



d d X + 



+ Q 



/S(W) 
/S(W) 



\ e 



+ 0,w) i {5 Q (W), J H }\ K + 
e.w) {^(w),j [q] }L 



The quantities 5{g[Q], J[ q ]}/5g k (6, W)|x; have the order 0(1) at e — > 0. We have also 

q %(^X),e^x(^X),...) ee ^M0',X), e cp x (0',X),...) 

according to relation (I3.16p . 

Let us introduce the functions 



^,., d )(<^x,...) = 



Using now relations (13. 28ft and f ]3.59j) we can write for the leading term of ( 13.631) 



%[Q]> J h) 



5g k (0,W) 



K[0] 



(-l) h+ - +l * k?(W)...k^(W) ... A£*(W)...*f'(W) 



d 



d 



o 



d 



h,—,ld 



de a i de a h de< qq< 



/S(W) 



Qi + O, WJ S^ ld) (cp(0, W), e<p w (0, W), . . . 



/C[0] 



?m(W) + 
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+ w«"(W) g M (W) Q M * + 0, (3.64) 



where 



^,., d) w), e^ w (0, w), . . = ^ + 0, s w , ^(w), . . ., fr^(wj) , . . .) 

Note one more property of the values 6{g[Q\, J{q]}/5g k (0, W)|^[o]- As we saw earlier, the values 
{<7[Q] 3 J[q]} are of the order 0(e) at e — >• on the submanifold K. This property is preserved also 
under the overall shift of the initial phase (I3.60p . 

Indeed, for 



¥>*(0,X) = $*^S + + A6>, S x , t/ 1 (X),...,?7 m+s (X)^ 



(3.65) 



we can write 



{#[Q]j ^[q]} = J Qi,6° 

x i {S a (Z),J [q] } $*^?l + + A0, S z , f/ 1 (Z),...,t/ m+s (Z)^) + 

+ / Qi + z ) C S"°) (* + + Sz ' ^( z )' • • • ' um+s ( z )) >•••)>< 

x ^(Z)(|^^ + 0(6) 

Because of the invariance under translations (13.601) the value of the bracket {^"(Z), J[ q ]} on the 
functions (13.651) is equal to its value on /C 

{S a (Z), J [q] } = e W -(Z)g„(Z) + 0(e 2 ) 
Similarly, we have on the functions ( 13.65)) 



U f^l + + A0, S z , U\Z), U m+s (Z)^j 



lo°"*(Z) $L ( ^ + + A0, S z , ^(Z), . . . , f/ m+s (z; 



We then obtain 



g,(Z)c^(Z) — [Q,(0,Z)^(0 + A0,Z)] + 0(e) = 0(e) 
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As a consequence, we can write 



f-2n 



2tt 



%[Q]> J [q]l 



5g k {0,W) 



<j) fe 



/Cfi 



SfWl 



0, S w , f/ 1 (W),...,f/ m+s (W)^) 



(27T) m 



d d W = 



(3.66) 

a = 1, . . . , m, for the main part of 5{<7[q], J^]} / 5 g k (6 , W) on /C. 

Using again the fact that relation (I3.66P contains arbitrary functions g At (W), appearing in the 
integrand expression in the form of local factors, we can rewrite (I3.66P in the stronger form 



2vr r2iT 



%[Q]' J [q]} 



"'J 5gk(0,W) 



<f> fc 



/C[0] 



SfW) 



0, S w , f/ 1 (W),...,f/ m+s (W)^) 



(27T)" i 



a = 1, . . . , m. 



At the same time we have 



${9[Q] > J[q\ } 



WW) 



0(e) 



%[Q]> J [q]} 



K 



0(e) 



K 



, vw 

(3.67) 



on the submanifold /C. 

We now turn back to the Jacobi identity (I3.57P for the functionals g\Q], J[q\, and Jy p y Using 
expansions of the values {J [q] , J [p] }, {g [Q] , J [q] }, {g [ci] , J [p] } in terms of # fc (0,W), S^W), f/ 7 (W) 
near /C, it is not difficult to see that after the restriction on K. the leading term (~ e) of relation 
( I3.57P can be written as 



y uqi^ (^ w )}i^] WW ) 



/C[l] 



(27r) m 



It n k (f) wUI ^Ig] ' J[q1 ^ 



(2 



7T 



- | {J M ,0*(0,W)} 



<H#[Q] , J [P}} 



m 5g k (e,w) 



d m e 

(27r) m 



tf'W = 



The above identity can again be written in a stronger form. Namely, making the change 
Qi(0,W) Q;(0, W)/i;(W), we get the corresponding change: Qi(0,W) Qi(0, W) ^(W), 
where /ij(W) are arbitrary smooth functions of W. From the form of the pairwise brackets of con- 
straints given by (13.391) . and (13.641) . it is easy to see then that the integrands are smooth functions 
of and W, containing /ij(W) in the form of local factors. By the arbitrariness of /Xj(W), we can 
omit the integration over W in the above integrals and write for every W: 
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K[l) 



d m e 

(27r) m 



+ 



/C[0] 



(2vr) m 



Finally, using relations ( 13 .401) and ( I3.67p . we can write the above identity as 



K[0] 



d m 9 



= 



d m 9 



27T /-27T 



JO 



A [i][p] I e 



+ 0,W 



8g k (0, W) 



/C[0] 



cTfl 
(2vr) m 



+ 



Sg k (0,W) 



K[0] 



d m 9 

(270™ 



where the functions ^^^(0, W) are introduced by formula (13.431) . 

Now assume that the values Af 1][q] (0, W) for U(W) G 5, according to (EHJ), can be represented 
in the form 



Af 1][q] (0,W) = -Bjg(W)B iW(1) (0,W) 

with some smooth in 6, 27T-periodic in each 6 a functions .Bj[ q ](i)(0, W). 
We can then write for U(W) G S: 



(3.68) 



JC[1] 



(2tt)' 



+ 



+ 



2?r /.2tt 



(! JO 

27T /.27T 



Jo 



[0][S] 



i[p](i) 



/S(W) 
/S(W) 



+ 0,W 



5^(0, W) 
<%[Q] > j [p]} 



d m 9 



(3.69) 



fy*(0, W) 



/C[Q] 



(27r) m 



where 
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h,...,l d x x 



/S(W) 



e, s w ,u 1 (w),...,u m+s (w) 

d d 



x A;" 1 (W)...A; 1 ' 1 (W) ... Q (W) . . . k d ' d {W) 



d 







d6 a i d6 a h 



Using expression ( I3.39P for the bracket of constraints on /C, as well as relations (I3.2ip . (j3.56p . 
f )3.62p . and the skew-symmetry of the operator BlL s ,(W), we obtain also the following relation: 



2tt f2ir 

ii Jo 



{g [Q] ,g (6,W)\\ m 6gk{e ^ w) 



d m e 



2tt i-2-k r 
ii JO 



/S(W) 



Ao]'[s]( W ) Qi [=— "• " 



S { J [q} > J [p]} 



d m 9 
jcpl] (2^F 



<f</*(0,W) 

Expanding also the remaining terms of identity (I3.69P according to (I3.64p . we get for U(W) G S: 



2n r 2ir 



JO 



^[s](W) Qj 



/S(W) 



0,W 



"H^q] > J [p]} 



5g k (e,W) 



oTO 



(3.70) 



27T /.27T 



d m 
o"'7o (2<r 



X 



X 



X 



£ kf(W) . . . k< (W) . . . kf(W) . . . k<(W) 

Odd 



■hf-,ld 

d 



d9 a h 



09 c 



Q l o.w) si" 



V e 



/S(W) 



K[0] 



+ 



+ 



2?T /■2-7T 



d m e 

o"'Jo Wr 



[o][sr 



\ e 



x 



X 



/S(W) 



il,...,i d 

9 



g 
3^ 



<9 



5 



00°? ' ' ' ^< 



Qi 



\ e 



e.w ^ w M0,w),... 



JC[0] 



P„(W) + 
= 
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Consider now the Jacobi identity of the form 

[9[P], [OIQ], Jin]}} + [9[Q], { J [q]^[P]}} + { J W> {9[P],9lQ]}} = (3.71) 

for arbitrary fixed functions q(X), and the functionals P(0,X) and Q(#,X) defined as before with 
the aid of arbitrary functions P(0,X), Q(#,X). 
According to the relations 

{<7 [Q , P] ,^(W)}L = 0(e) , {J [q] ,<7 fe (0,W)}L = 0(e) , { J [q] , U»(W)}L = 0{e) 



it's not difficult to see that after the restriction on K, the major term (in e) of (13.711) will be written 

as 



In n k (f) WUI 5 ^ 9lQ] I J[q1 ^ 

{g [P] ,g (0,W)\\ m Snk(Q 



d m 9 
(2vr) m 



d d W - 



Sa n k (fi WUI ! ^ 

{g [Q] ,g (0,W)\\ m Sgk{6 ^ w) 



d m e 

(2tt)™ 



<W = 



Again recalling that qv(W) are arbitrary functions of W appearing in the integrand in the form 
of local factors, we can write the above relation in a stronger form. That is, for every W 



2vr i-2-k 



JO 
2tt 1-2-e 



in WUI 6 ^ 9[Q] ' ^ < '"' f> 

{g [P] ,g (0,W)\\ m 6gk(0 ^ w) 



Sn n k (f) WUI ^[P] ; J [q]} 



m (2t) w 

d m e 



(3.72) 







As well as in the case of identity (13. 70 p . using relations (13. 39 p . (I3.2ip . (I3.56p . (I3.67p . we can write 
identity (I3.72p in the form: 



2n ,2tt d m 6 

"'7 (2<P 



^[S](W)P,f^-: O.W 



(3.73) 



x 



x 



(-1)^+-+^ kf 1 (W)...kf 1 (W) ... kf{W)...kf d {W) x 

h,—,ld 

d d d d 



06< d9 a h d$ a i 90 a t 



(h(^ + e,w)s^ ld) (cp(e,w),...) 



V e 



K[0] 



?m(W) + 



2tt /.2tt 



^0 



d m 
(2vr) m 



\ e 
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X 



X 



.'l id 

d 



H hid 



fc° 1 (W)...A; 1 !l (W) ... ^ 1 (W)...^' d (W) x 















d6 a i QQ a h 



d6 c 



de c 



Pi 



/S(W) 



o^)s% iUd) (^,W),...) 

/S(W) 



+ w^(W) g M (W) P, 



V e 



f 0,W 



?m(W) + 
= 



Note now that the values of Q(0, X) and P(0, X) are arbitrary 2-7r-periodic functions of d, satis- 
fying conditions (13.561) . In particular, we can put in (I3.73P at U(W) G S 



P(0,W) = B [p](1) (0,W) or P(0,W) = B [q](1) (0,W) 
By analogy with (I3.64p we introduce for convenience the notation for U(W) G S: 



<H0[eB [p](1) ] > J [q]l 



<5# fe (0,W) 



(3.74) 



(3.75) 



£[1] 



£ (-1)^-^ ^(W)...^(W) ... fc^(W)...fc^(W) 



<9 



9 



<9 



9 



de c 



BiMH) (^ + ^w)^ i .. id) (^,W) 
+ W a "(W)?,(W)B fcIp ] (1) , s , 



S(W) 



<k(W) + 



e.w 



for arbitrary smooth functions q(X), p(X). 

Note that the functional <?[ e B [p](1) ] is not defined on the whole functional space, so relation (13.751) 
plays just a role of a formal notation for U(W) G S. 

Using now (I3.73P for the functions (13.741) . we can rewrite (13.701) in the form 



2vr p2n 



d m 9 

'"Jo Wr 



B^ ][s] m Qj 



S(W) 



0,W 



x 



x 



\ Sg k (6,W) 



S{g[e-B lp](1) ] , J[q]} 



K[l] 



Sg k (0,W) 



+ 



%[eB [q](1) ] , J[p]} 



K[l] 



6g k (0,W) 



provided that U(W) G S. 



Using the skew-symmetry of the Hamiltonian operator B^ S ^(W) on A we can then write 



■2jt 



" d "' e Q] (^l + e,w 



X 



n (27T 
B jk 



<K J [q] , ^[p]} 



K' S]( W ) \ 6g k (0,W) 



S{g[eB [p]w ] , J[d]} 



K[l] 



5g k (e,W) 



<H#[eB [q](1) ] , J[p)} 



K[l] 



5g k (0,W) 



= o 
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The values Qj(0, W) are arbitrary smooth 27r-periodic functions of with the only restriction 
(I3.56p . We know also that the value in the brackets is a smooth 27r-periodic in each 8 a function of 
for U(W) G S. As a consequence, we can write for U(W) G S 



BiLm 1 HJ, " , - J, ' ,] 



5g k (0,W) 



${9[eB lp]m ], J[q\} 



/C[l] 



Sg k (0,W) 



H9[eB mi) ],J\p]} 



K[l) 



8g k (0,W) 



X>f q>p] (U(W),U w (W)) % a 



a=l 



S(W) 



0,U(W) 



with some coefficients a[* qp ](U(W),Uw(W)). 

The values in parentheses are smooth 27r-periodic functions of at U(W) G S. According 
to Lemma 3.6 we can therefore say that up to a linear combination of the regular eigenvectors 
v [u(w)](S(W)/ e + 0) of J3jgf [sl (W) = i^ sl (U(W)), corresponding to zero eigenvalues, the value 



in parentheses is a linear combination of the variation derivatives f l2.10p . generating linear shifts of 
the phases on A. For a complete Hamiltonian set of the integrals (J 1 , . . . ,1) we can then write 
according to fl2TTTj) and f!3"3Ij) 



<K</[ q ], J[ P ]} 



5g k (0,W) 



%[eB [p](1) ], ^[q]} 



K[l] 



5g k (0,W) 



H9[c-B M{1) ], J\p]} 



K[l] 



m+s 



= EWi(u(w),u w (w))4«) I 



(W)] 



9=1 



0,W 



S(W) 



/C[l] 







(3.76) 



with some coefficients 6 9[q>p] (U(W), U W (W)) at U(W) G S. 
Consider now the Jacobi identity of the form 

{{J[q\, J[p]}, J[r)} + {{ J [ph J [r]} , J[q\} + { { J[r] , J[q\ } , J[p] } = 

with arbitrary smooth functions q(X), p(X), r(X). 

In the main (~ e 2 ) order on K. the given identity leads to the relations 



H j Wi J {p]} 



+ 



5S a (W 

5{J[q], J[p]} 



K[l] 



{S a (W),J [r] }\ K[1] d d W + c.p. + 



+ 



SU-r(W) 

H J {dh J[p]} 

Sg k (0,W) 



/C[l] 



{U^(W),J [r] }\ m d d W + c.p. + 



(3.77) 



j igW ,j lr]} \ m — 



d d W + c.p. 



(a = 1, . . . , m, 7 = 1, . . . , m + s). 

Again, using relations (I3.40p . (I3.43j) and (I3.62p . we can replace identity (I3.77P by the following 
relation 



H J [<ti> J[p]} 



5S a (W) 



K[l] 



{S a (W), J [r] }\ K[1] d d W + c.p. + 
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+ 



H J lib J [p]} 



+ 



5U-r(W) 
2W f 2n H J [nh J [p]} 

---J Sg k (e,W) 



K[l] 



{U\W),J [r] }\ K[1] d d W + c.p. + 



A 



K[l] 



MM ^ 



/S(W) 



0, W 



(3.78) 



d m 9 
(2vr) m 



d d W + c.p. = 



Using relations (I3.44p and the representations f!3.68j) and (I3.76p . we can write for U(W) G S: 
J "'J 5gk(6,W) 

d' ,! () 



Ak /S(W) 



d m 

+ e,w\ + c.p. 



(3.79) 



2?r /-27T 




'0 JO 

»2tt p2ir 



/S(W) 



^(W) -0.W 



%[ e B [q](1) ], J[p]} 




+ /... 

'0 ^0 



^ [ s](W)fi j[q](1 )(^ + ^W 



^ fc (^ ; w ) 

%[eB [p](1) ], </[«•]} 



<H#[ e B [p](1) ], J[q]} 



n2?r r 



^ [S] (W)S i[p](1) (^ + 0,W 



5g k (6,W) 

%[eB [r](1) ], J[q]} 



^ WW) 

_ <H#[eB [r](1) ], </[p]} 

/C[l] 



6g k (6,W) 



W,W) 

%[eB [q](1) ],^[r]} 



/C[l] 



/C[l) 



JC[1] 



<fc/*(0,W) 



(2?r) m 
(2vr) m " 



AC[1] 



(27T) 



Similarly to earlier arguments, substituting now in identity ( 13. 73 ft : 

Q(0,W) = B [r](1) (0,W) , P(0,W) = B [p](1) (0,W) 
for U(W) G 5 we obtain the identities 



2n r2n 



B [o][s]( W ) B J[P](D 



^ S(W) 



+ 0,w 



5{c/[ £ B [r](1) ],J[q]} 



V=(0,W) 



d m e 



(3.80) 



2tt /.2tt 



^ [S] (W) B j[m ^l + e,w 



%[eB [p](1) ], J[q]} 



fy fc (0,W) 



d m e 



o 



Using the cyclic permutations of the functions q(X), p(X), and r(X) in identity ( I3.80p . it's not 
difficult to see that the right-hand part of relation (I3.79j) is identically equal to zero at U(W) G S. 
It's not difficult to see also that the left-hand side of the relation (13.791) is a smooth regular function 
of the parameters U(W) = (Sw, f/ 1 (W), . . . , U m+s (W)) and their derivatives. Using the fact that 
the set S is everywhere dense in the parameter space U, we can conclude that the left-hand side of 
equation (13. 79j) is identically equal to zero under the conditions of the theorem. 

We have, therefore, that under the conditions of the theorem, the identity (13.771) implies the 
relations 



H J [q\i J \p]} 



+ 



Using the relations 



5S a (W) 

<H J [q]> J [p]} 



K[l] 



{S a (W),J [r] }\ K[1] d d W + c.p. + 



K[l] 



{tr(W), J [r] }\ K[1] d d W + c.p. 
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H J l<ib J [p}} 



6S a (W) 



H U lq\> U \p] 



DN 



K[l] 



5S a (W) 

{trr(w),j w }| 



Hu M ,u [p] } 



DN 



K[l] 



{Ui(W) 7 U [r] } 



DN 



we get then the identities 



+ 1 ^;^ {tr > (W ),g M } IM , rty + c.p. S (3.81) 

(summation in a = 1, . . . , m, 7 = 1, . . . , m + s). 

Here the notation {^"(W), Um}dn means the pairing of S" a (W) as the functional of U(Z) with 
U[ r ] given by the Dubrovin - Novikov skew-symmetric form and coinciding with the corresponding 
value for the bracket ( 13.551) . Easy to see then that relations (13.811) give in particular the Jacobi 
identities for the bracket fl335l) on the space of fields (S Q (X), U\X.), . . . , U m+s (X.)). 

Theorem 3.1 is proved. 

We will not consider here applications of Theorem 3.1 in details. The example of the application 
of Theorem 3.1 in one-dimensional case can be found in [25]. In general, we can say that Theorem 
3.1 gives justification of the procedure of the bracket averaging for a wide class of systems having 
multi-phase solutions. 

Using Lemma 3.5 we can formulate now the theorem that gives the justification of the procedure 
of the bracket averaging for a regular Hamiltonian family of single-phase solutions of system (12.51) . 

Theorem 3.1'. 

Let A be a regular Hamiltonian family of single-phase solutions of system Ii2.5\) . Let (J 1 , . . . , I ) 
be a complete Hamiltonian set of commuting first integrals of system Ii2.5\) on A having the form 
Ii2.8\) . Then the bracket 

{S(X),S(Y)} = , {S(X) , ET(Y)} = u> (S x , ^(X), . . . , U S+1 (X)) <f(X-Y), 
{WpQtU'OO} = (AJLo) (Sx, ^(X),...,^ 1 ^)) <mx-y) + ...+ 
+ (AZ. 01 )(Sx,U 1 (X),...,U s+1 (X.)) 5 xd (X-Y) + 

+ [(Q 7PP >(Sx, U 1 (X),...,U S+1 (X))] XP 5(X-Y) , 7> p=l,..., a + l (3.82) 

on the space of fields (S(X) , f/ 7 (X)) ; 7 = 1, . . . , s + 1, satisfies the Jacobi identity. 

Let us prove now the second theorem justifying the invariance of the bracket averaging procedure. 
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Theorem 3.2. 

Let A be a regular Hamiltonian family of m-phase solutions of A2.5\) . Let (I 1 , . . . , I N ) and 
(J 71 , . . . , I' N ) be two different complete Hamiltonian sets of commuting first integrals of A2.5\) hav- 
ing the form $2.81) . Then the brackets $3. 55\) obtained using the sets (J 1 , . . . , I N ) and (I' 1 , . . . , I' N ) 
coincide with each other. 

Proof. 

Note that the sets (I 1 ,. . . , I N ), (I' 1 , . . . , I' N ) correspond to two different systems of coordinates 
(U 1 , . . . , U ), (U 11 , . . . , U' N ) on the family A, given by the averages of the functionals I and I'. Let 
us first prove that the Dubrovin - Novikov form obtained using the set (J' 1 , . . . , I' N ) coincides with 
the form, obtained using the set (J 1 , . . . , I ), after the corresponding change of coordinates. 

U' v = U' V (U) 

Consider the functionals 

t>2ir p2ir d m 9 

, u = l,...,N 



PZ7T P Z7T 

J"(X) = ... P»(<p,ecp^... 
Jo Jo 



(2vr) m 

Consider now the values of the functionals S a (X), 7 7 7 (X), and the constraints g l (6, X), introduced 
with the aid of the functionals J(X), as a coordinate system in the neighborhood of the submanifold 
K. 

We obviously have the relations 

r{X)\ K = W(X) + 0(e) , r(x)\ K = E7*(X) + O(e) 
on the submanifold /C. 

For values of the functionals J' V {X) on the submanifold K, we can then write the relations 

r{x)\ K = u'»(j(x)) + ^j^(B Xl u 1 (x),...,u m ^(X) 1 ...) 

l>l 

where are smooth functions of (Sx, f^ 1 (X), . . . , U m+s (X), . . . ) and their derivatives, polynomial 
in the derivatives (Sxx, Ux, • • • ) and having degree I according to our previous definition. 
Expanding the values of J' U (X) in the neighborhood of the submanifold /C, we can write 

r(X) = U'»(J(X)) + ^e i jS(Sx,?/ 1 (X),... l tr n +'(X),...) + 

l>\ 

*2n r2n d mg 



+ //.../ lf(X,0,Y,e) <f(0,Y) 
J Jo Jo 



d d Y + 0<~ 2 " 



(27T) 

where, according to the form of constraints (13. lip , we can put 

If(X,0,Y,e) = Yl nr (/l - W (^,X), e ^x(^,X),...)^ e ' 1+ -^ S hX ^ ldXd (X-Y 



h,—,ld 



E J& M^X),e^ x (^X),... 

h,...,l d U ^hx\..l d xd 



e tl+ - +td <W...^(X-Y) 

K 
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Considering the functional J' q j = f q v {X) J' U {X) d d X with arbitrary smooth (compactly sup- 
ported) functions q u (X), we can write in the vicinity of /C: 



JL 



?,(x) 



U w (J(X)) + e ' J(0 ( S *' ^(X), • • • , U m+s (X), ...) 



1>1 



d d X + (3.83) 



h+-+ld e h+-+ld 



h,...,ld 



d h d ld 



T Qu{Y) n;^- W (ip(0,Y),eip Y (0,Y),... 



dYl 1 dY c 



X 



d m e 



x 9\0,Y) j^-J'Y + 0(g 2 



The leading term (in e) in the second part of expression (13.831) is given by the expression 



QuOO £ (-l) h+ - +ld k^(Y)...k^(Y) ... k a /(Y)...k d l *(Y) x 

h,...,ld 

x n«*-w , „ (• (^ + », y) , . . . ) 9 < (9, Y) ™ 
h...e a i...e id V V e J J ( 27r ) 



and coincides with the value 

"2tt /.2tt 



JO 



where 



q,(Y)C 



Jo Jo 



i[U(Y)] I e 



d m 

+ e)g\e,Y) -^—d'Y 

(2tt) 



(2vr) m 



<p(0)=*(0,U) 



As the values Qui(^)' ^ ne va l ues Qui(^) represent regular left eigenvectors of the operator L 1 -^ 
corresponding to the zero eigenvalue. In the case of a complete regular family of m-phase solutions, 
we have therefore 



for some functions r^(U). 

We can write, therefore, up to quadratic terms in g(0, X): 



J q v {X)J' v {X)d d X = J q v \ 



X) 



E/*(J(X)) + J2e l j'^{S^U\X),...,U m+s (X),...) 



i>i 



d d X + 



+ / g„(x) 



* + 0(6) 



^' X) + ° (g2) 
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Consider the Poisson brackets: 



Wd'-WI/c = {/^(X)JnX)rf d X, j p„{Y)J>»{Y)d d Y 



K 



dU w 
^ (X) 9^ (X) 



J A (X) , | p„(Y) J'^Y)d d Y 



a: 



+ |g,(X)^rj(U(X)) 



(<z) 

i[U(X) 



S(X) 







0(e) 



By Lemma 3.2' we have the relation {g l (6, X) , ■/r p i}|x: — 0(e) on the submanifold K. In addition, 
completely analogous to relation (I3.44p holds the relation 



2tt r 2n 



JO 



K 



(?) 

i[U(X)] 



SfX) 



0j (<f(0,X), Jf p] }k[i] 



eP0 
(27r) m 







by virtue of the original dependence of the constraints g % (0, X). We thus obtain: 



J j q v (X) {r(X),J'»(Y)}\ K p,(Y)d d Xd d Y 



d d X + 0(e 2 ) 



= / fc,(X) ^(X) { J A (X) , | P/i (Y) J'^(Y) d*r} 
Repeating the arguments for the functional j p^{Y) J'^iY) d d Y we finally obtain 

j j q u (X) {<T(X), J'»(Y)}\ kPii (Y) d d Xd d Y = 



(3.84) 



dU 



ff „(X) ^(X) (j A (X), r(Y)}\ K _(Y)p M (Y) d d XcfV + 0(e 



<9£/ c 



Given that the principal (in e) terms in the expressions { J A (X) , J a (Y)}\ic and { J /U (X) , J' m (Y)}|a; 
coincide with the Dubrovin - Novikov forms, obtained with the aid of the sets (I 1 , . . . , I N ) and 
(I' 1 , . . . , I' N ) respectively, we conclude from (13.841) : 



{^(x),^(y)>; 



DN 



(X){ f / A (X),^(Y)} r .^-(Y) 



dU^ v 11 v ,iDN dU a 

which means the coinciding of the forms {...,... } DN and {...,... }' DN . 
As a result, we can claim that the relations 



{A£(U(X)) , *f(U(Y))} = , {^(U(X)),fT(Y)} = [^(U(X)) 5(X — Y)] 
(a, /3 = 1, . . . , m, 7 = 1, . . . , m + s) and the expressions for 



XP 
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{IP(X),U"(Y)} DN , j,p=l,...,m + s 
transform into the expressions 

{^(U'(X)),^(U'(Y))} = , {^(U'(X)),[^(Y)} = [u/^(U'(X)) 5(X-Y)] XP 

and 

{Z7*(X) , U">(Y)}' DN 

after the change of coordinates 

(U\...,U N ) (U'\...,U' N ) 

on A. We can claim then that brackets (I3.55P obtained with the aid of the sets (J 1 , . . . , / ) and 
(J' 1 , . . . , I' N ) transform into each other after the change of coordinates 



(^(X), . . . , S m {X), U\X), U m+s {X)) -> (S\X), . . . , S m {X), U a {X), U' m+s {X)) 
where U'\X) = W (S x , U^X),..., U m+s (X)). 

Theorem 3.2 is proved. 

Let us note that Theorem 3.2 means, in particular, that bracket f!3.55|) is also invariant with 
respect to the choice of the functionals (J 1 , . . . , J m+S ) among the full set {I u , u = 1, . . . , N}. 

Finally, we prove the theorem about the Hamiltonian properties of the Whitham system (I2.40p - 
(I2.4ip under the same conditions as before. 

Theorem 3.3. 

Let A be a regular Hamiltonian family ofm-phase solutions of A2.5\) . Let (J 1 , . . . , I N ) be a complete 
Hamiltonian set of commuting first integrals of 112. 5\) having the form A2.8\) and H be the Hamiltonian 
function for system A2.5\) having the form [2. Sty) . Then the Whitham system \2.4ty - j2.41\) 
Hamiltonian with respect to the corresponding bracket $3.55\) with the Hamiltonian function 



IS 



H av = J (P H ) (S x , ^(X), . . . , U m+s (X)) d d X 



Proof. 

By Theorem 3.2, without loss of generality we can assume that the Hamiltonian functional H be- 
longs to the set (J 1 , . . . , / m + s ) ; H = J Mo . It is easy to verify then that the corresponding Hamiltonian 
H av generates in this case the system 

S« = uj a ^ (S x , U 1 ,..., U m+S ) , a = 1, . . . , m 

W T = (Q^xi + ... + (Q^ od ) xd , 7 = l,...,m + s 
i.e. exactly system fl2^0|) - (12TU1) . 

Theorem 3.3 is proved. 
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4 On the canonical forms of the averaged brackets. 

Let us consider now " canonical forms" acquired by the averaged brackets in some special coordinates. 
We first prove here the following lemma: 

Lemma 4.1. 

Consider the averaged bracket fP5)) in the coordinates (S^X), . . . , S m (X), E/^X), . . . , t/ m+s (X)). 
There exists locally an invertible change of coordinates 

(S 1 (X), S m {X), U\X), U m+s {X)) ->• 

-> (^(X), . . . , S m (X), Q 1 (X), . . . , Q m (X), N\X), . . . , N l (X)) 

where 

Q a = Q a {Sx,U\...,U m+s ) , N l = N l (Sx,U\...,U m+s ) 
such that the bracket A3. 55\) has in the new coordinates the form: 

{S a (X),S^(Y)} = , {S a (X),Q p (Y)} = ^(J(X-Y) , {S a (X),N l (Y)} = (4.1) 



{Q a (X),Q p (Y)} = J a/3 [S,N](X-Y) , {Q a (X),N l (Y)} = J< [S, N] (X - Y) 

{N l (X) , N q {Y)} = j' 9 [S, N] (X - Y) (4.2) 

(a, (3 — 1, . . . , m, I, q — 1, . . . , s), where J a p, J l a , J lq are local distributions of the gradation degree 1 
given by linear combinations of the function S(X — Y) and its first derivatives with local coefficients, 
depending on the values (Sx, N(X), Sxx, Nx(X)). 

Proof. 

Let us fix the parameters k® = S^p on the family A and consider the coordinates 
U = (U 1 , U m+S ) on the (m + s)-dimensional submanifolds k® = const (for all a,p). Let us 
consider the vector fields 

§«■) = (u; al (Sx, U),...,u am+s (Sx, U)Y 

on the submanifolds k^ = const. 
From the Jacobi identities 

{{^(X),S Q (Y)},^(Z)} - {{U v (X),S p (Z)} ,S a (Y)} = 
for the bracket (13.551) it is easy to get the relations 



f («) ' 



a, (3 = 1, . . . , m 



for the commutators of the vectors fields £( a ) on the submanifolds = const. 
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According to relations ( 13 .2p for the frequencies generated by the functionals (J 1 , . . . , J m+S ) we 
can state also that the set {£(«)} is linearly independent at every point. 
We can claim then that on every submanifold 

{kp = const , a — 1, . . . , m , p — 1, . . . , dj 
there exists locally an invertible change of coordinates 

(U\...,U m+s ) -> (Q 1 (S x ,U),...,Q m (S x ,U), iV 1 (S x ,U),...,iV s (Sx,U)) 

depending smoothly on the parameters = Sxp, such that the coordinate representation of the 
vector fields £( a ) on these submanifolds has the form 

= (1,0,...,0)*, ... , gm) = (o,... ,0,1,0,. ..,0)* 
It is not difficult to see then that the corresponding change of coordinates 

(S(X),U(X)) (S(X),Q(X),N(X)) 

provides relations (14. ip for bracket (I3.55p . 

The general form of the distributions J a p, J a , J lq follows then from the form of bracket ( I3.55P 
and for the proof of the Lemma we have just to prove the relations 

5 J a p 1 5Q 7 (Z) = , 5J l a /5Q y (Z) = , 5J lq / SQ 7 {Z) = , 7 = l,...,m (4.3) 

for the functionals J a p, J l a , J lq . 
Using the Jacobi identities 

{{Q a (X),Q^(Y)},5T(Z)} + c.p. = , {{Q a (X) , N l (Y)} , £T>(Z)} + c.p. = 

{{N l (X),N\Y)} ,^(Z)} + c.p. = 
we then easily get relations (14. 3p . 

Lemma 4.1 is proved. 

Let us note here that according to their definition the coordinates (Q a , N l ) are defined modulo 
the transformations 

Q a -» Qa + g«(S x , N) , N l -)> iV z (S x , N) 
where the transformation N l — > N l (Sx, N) is invertible for every fixed Sx- 

Let us consider now a special case, when the additional parameters (n 1 , . . . , n s ) are absent on the 
family A and the full regular family of m-phase solutions of system (12. 5p can be parametrized by the 
mid + 1) parameters (hp*, u a ), a = 1, . . . , m, p = 1, . . . , d. 

In this situation we can prove the following theorem about the canonical form of the averaged 
bracket f l335|) . 



57 



Theorem 4.1. 

Let system (12. 5\) be a local Hamiltonian system generated by the functional A2.33\) in the local field- 
theoretic Hamiltonian structure 112. 3 2\) . Let A be a regular Hamiltonian family of m-phase solutions 
of (12. 3\) and (J 1 , . . . , J m ( d + 1 )) ft e a complete Hamiltonian set of commuting integrals \2. 8\) for this 
family. 

Let relations (\3. 55\) represent a Poisson bracket on the space of fields 

(^(X), . . . , S m (X), U\X), U m (X)) 
Then locally there exists an invertible change of coordinates 

{S\X), S m (X), U\X), U m (X)) (^(X), . . . , S m (X), Q^X), Q ra (X)) 

where Q a = Qa(Sx, U) such that the bracket (13. 55\) acquires in the coordinates (S(X) , Q(X)) the 
following non- degenerate canonical form: 

{S a (X),S?(Y)} = , {S a (X),Q p (Y)} = 5%5(X-Y) , {Q a (X) , Q P (Y)} = 
Proof. 

Let us give here just a proof in abstract field-theoretical form. The same considerations can be 
made also in the differential-geometric form after concrete substitutions for the coordinate transfor- 
mations. However, we would like to skip here the corresponding calculations to avoid rather bulky 
expressions in the text. 

In accordance with Lemma 4.1 let us write down the averaged bracket f!3.55j) in coordinates 
(S a (X) , Q«(X)) such that 

{S a (X),S^(Y)} = , {^(X^Q^Y)} = 6%6(X-Y) 

{Q Q (X),Q /3 (Y)} = J Q/3 [S] (X - Y) 

From the Jacobi identities 

{{Qa(X),Q/3(Y)} ,Q 7 (Z)} + c.p. = 
we then easily get the relations 

5J aP [S](X,Y) ^ 7 [S](Y,Z) 5J ja [S}(Z,X) = 

5S-y(Z) 5S a (X) 5SP(Y) ~ 1 ' ' 

for the form J afj [S](X, Y) = J Q/9 [S](X - Y). 

Relations (14 Ah mean that the 2-form j aj g[S](X, Y) on the space of fields (S' 1 (X), . . . , S m (X)) is 
closed and can be locally represented as 

J rql , Y v ^ Sg a [S)(X) _ %[S](Y) 
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for some 1-form g a [S](X) on the same phase space. Easy to see now that the change 

Qa(X) = Q a {X) - g a [S](X) 

gives the necessary coordinates Q a (X). 

From the form of the functionals J aj g[S](X, Y) it's not difficult to show also that the 1-form 
g [S](X) can be chosen in the form g a [S](X) = g a (Sx) with some smooth functions g a (Sx), which 
gives the necessary form of the corresponding coordinate transformation. 

Theorem 4.1 is proved. 

Let us note that in more general case of the presence of additional parameters (n 1 , . . . ,n s ) we 
can expect in fact, that the bracket (14. ip - (14. 2p can have more complicated, maybe even degenerate, 
form. 

We can note also that the bracket {iV'(X) , N q (Y)} given by (14. 2p represents a Poisson bracket 
on the space of fields N(X) for any fixed values of S(X) as follows from the form of the bracket (14. ip 
- (14. 2p . This actually gives rather strong restrictions on the form of the bracket 



At last we consider just a very simple example of the averaging of a non-degenerate bracket in 
the single-phase case. Let us consider the nonlinear wave equation in d spatial dimensions having 
the form 

<p tt - A<p = V\ip) (4.5) 

with some potential function V(ip). 

The periodic one-phase solutions of (14. 5 p are defined by the equation 

_ k 2 fc 2) ^ = y l{jp) (4 

which implies 

^— A_ 54 ^ = V{(p) + C 

with some integration constant C. We can assume for example that we have here u 2 > k 2 though 
it does not affect the final conclusions. 

The solutions if (9, u, hi, . . . , kd) are then given by the quadrature 



2 J v>V(tp)+C 

where the value of <p oscillates between two subsequent zeros of the expression V(<p) + C, while the 
condition 



" 2 - fc i ^ I - j£ = 2n 

2 J y/V(ip)+C 

fixes the value of the integration constant C. 

Let us denote the corresponding function (p(9, u, ki, . . . , kd) as $(# + 9 , u, ki, . . . , kd) 
Equation (14.51) can be written in the Hamiltonian form 
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with the Hamiltonian function 



H = I Puis, r....)d d x = 
and the Poisson bracket 



i> t = Aip + V'(<p) 

r (Vy) 2 
2 2 



- V(<p) d d x 



(4.7) 



{<p(x) , <p{y)} = , {^(x) , V(y)} = 5(x - y) , {^(x) , V(y)} = o 
The family A is represented then by the family of the functions 



<P 



$(# + #o, uj,ki,...,k d ) 
u$o(9 + 9 , u,ki,...,k d ) 



parametrized by the values (u,k\, . . . , ka) and the initial phase ho- 
using equation (14. 6 p it is not difficult to show that the family A is a complete regular family of 
single-phase solutions of ( 14. 7ft . The set of the commuting functionals (J 1 , . . . , I d+1 ) is given here by 
the momentum functionals 



P q {ip, ^,...)d d x 



ip x g ijj d d x , q — 1, . . . , d 



and the Hamiltonian functional (14.81) , which give a complete Hamiltonian set of the first integrals for 
the family A. The values of the averaged densities (P q ) on A are equal to 



(P q ) 



1^*2 dd 
o 9 2^ 



kg Q 



where 



Q 



2tt 



de 

2^ 



while the functionals I q generate on A linear shifts of 8q with the frequencies u q = k q . 
The pairwise Poisson brackets of the densities P 9 (x), Pff(x) can be written in the form 

{P,(x),P(y)} = P 9 (x)Mx-y) + Pz(x)4,(x-y) + [P g (x)] sI «J(x-y) (4.9) 



{P 9 (x) , P H (y)} = ^ 2 (x) ^(x - y) + <p xq ^ M x _ y) + 

i 

5(x-y) + ^ iVxiV x i) x i 5(x-y) 



+ 



~ \ (V^) 2 + V{<p) 



(4.10) 



{p„(x) , Pn(y)} = 2 - y) + E [ p <( x )L< J ( x - y) (4- 11 ) 

i I 

Thus, if we choose the coordinates 5(X), [/( 9 )(X) = (P 9 )(X) for some g = the corre- 

sponding bracket (13.821) will be written in the form: 
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{S(X) , S(Y)} = , {S(X) , U (q) (Y)} = S xq S(X - Y) 

(4.12) 

{U {q) (X) , U {q) (Y)} = 2U {q) (X)5 xq (X-Y) + U {q)x «5(X-Y) 
It's not difficult to check that in the coordinates 

S(X) , Q(X) = (P 1 )/3 X1 = ... = (P d )/S xd 
any of the brackets (I4.12p takes the form 

{S(X),S(Y)} = , {S(X),Q(Y)} = 5(X-Y) , {Q(X) , Q(Y)} = (4.13) 

We can see then that all the brackets (I4.12p represent in fact the same bracket (I4.13P in different 
coordinates. Thus, we can claim that all the brackets (I4.12p transform into each other after the 
corresponding change of coordinates U( q )(X) = S Xq Un\(X)/S x i. 

The averaged density (Ph) on A can be represented in the form 



(Ph) 



2tt 



co 2 + k\ + • • • + k 2 a 



2- 



% 27T 



d6 
2^ 



2t 



■~k 



dB_ 
2^ 



Using the fact that the variation derivative of the second integral is equal to zero on A according 
to (14.61) it is not difficult to obtain the relation 



k, 



d(P H ) = udQ + QS^-dh 



(4.14) 



From (I4.14p it is easy to get the following relations in the bracket (14 . 1 3 [) : 

{S{X) , {P H )(Y)} = uj(S x ,Q)5(X-Y) 

{(P H )(X) , (P H )(Y)} = 2^<PA(X)MX-Y) + J2( p i)* 5 ( X "Y) 

i i 

From expression (14.111) for the pairwise Poisson brackets of the functionals P^(x), Pjf(y) we can 
see then that the bracket (I3.82p obtained with the aid of the functional H coincides with the bracket 
(I4.13P (and all (14.121) ) after the corresponding change of coordinates. 

The Whitham equations can be easily written using the averaged Poisson bracket and the func- 
tional H av = J (Ph) d d X . Thus, in the coordinates (S(X), Q(X)) the Whitham system takes the 
form 



w (Sx, Q) 



Q 



T 



E 



S x i 

UJ 



Q 
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